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FOREWORD 


This report describes a dynamic pressure data base and data 
base management system developed to characterize the Space 
Shuttle Main Engine (SSME) dynamic pressure environment. 
The data base represents dynamic pressure measurements 
obtained during single engine hot firing tests of the SSME. 
Software is provided to permit statistical evaluation of 
selected measurements under specified operating conditions. 
An interpolation scheme is also included to estimate 
spectral trends with SSME power level. This report was 
prepared by Wyle Laboratories Scientific Services and 
Systems Group for the National Aeronautics and Space 
Administration, George C. Marshall Space Flight Center. The 
work was performed under NASA contract NAS 8-34343, entitled 
"Flow Dynamic Environments in High Performance Rocket 
Engines . " 

The author wishes to acknowledge the contribution to this 
study by Dr. G. Meares, Chief Architect of the data base 
management software described herein. Messrs. B. Dobbs and 
D. Duck contributed long hours over a digitizer pad to 
accomplish data base input. The singular value 
decomposition software for spectrum interpolation was 
developed by Dr. J. Jong. Mr. T. Nesman, MSFC technical 
contract monitor, provided continuing support through 
informal project reviews and served as a focal point for 
definition of SSME data requirements. 
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TIi m i ypnnrt — .-itr ■nn fees- a_ dynamic pressure data 
baseQmanagement system for measurements obtained from space 
shuttle main engine (SSME) hot firing testsAsfe '* 
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time histories, and power spectra of the dynamic pressure 
measurements at selected times during each test.- Test 
measurements and engine locations are defined along with a 
discussion of data acquisition and reduction procedures. A 
description of the data base managemenl^analysis system is 
provided and subroutines developed for obtaining selected 
measurement means, variances, ranges and other statistics of 
interest are discussed. A summary of pressure spectra 
obtained at SSME rated power level is provided for 
reference. Application of the singular value decomposition 
technique to spectrum interpolation is discussed and 
isoplots of interpolated spectra are presented to indicate 
measurement trends with engine power level. Program 
listings of the data base management and spectrum 
interpolation software are given. Appendices are included 
to document all data base measurements. 


ABSTRACT 


This report describes a dynamic pressure data 
base/management system for measurements obtained from space 
shuttle main engine (SSME) hot firing tests at the National 
Space Technology Laboratories and the Santa Susana test 
facility operated by Rocketdyne. The data were provided by 
NASA/MS FC in terms of engine power level and r.m.s. pressure 
time histories, and power spectra of the dynamic pressure 
measurements at selected times during each test. Test 
measurements and engine locations are defined along with a 
discussion of data acquisition and reduction procedures. A 
description of the data base management/analysis system is 
provided and subroutines developed for obtaining selected 
measurement means, variances, ranges and other statistics of 
interest are discussed. A summary of pressure spectra 
obtained at SSME rated power level is provided for 
reference. Application of the singular value decomposition 
technique to spectrum interpolation is discussed and 
isoplots of interpolated spectra are presented to indicate 
measurement trends with engine power level. Program 
listings of the data base management and spectrum 
interpolation software are given. Appendices are included 
to document all data base measurements. 
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SECTION I 

INTRODUCTORY SUMMARY 


The Space Shuttle Main Engines (SSMEs) are designed to 
provide maximum performance within stringent constraints on 
size, weight, and efficiency. As a result the SSMEs are 
required to operate under extreme temperatures, with high 
fluid pressures and rotational pump speeds. Since the Space 
Shuttle Vehicle (SSV) is man-rated, and a key to the cost 
effectiveness of the SSV concept is hardware reusability, 
system reliability is of paramount importance. 

Even under nominal operating conditions, the rocket engine 
systems and components are subject to significant pressure 
fluctuations. These pressure fluctuations emanate from 
pumps, turbines, valves, and in some instances highly 
turbulent propellant line flows at high dynamic pressures. 
These dynamic environments have temporal characteristics 
such that in some cases they appear as complex periodic 
functions and in other instances as wide-band random noise. 
Intense pressure oscillations have been a source of high 
cycle fatigue damage and SSME component failure. 
Measurement and empirical classification of the dynamic 
pressure environment is therefore desirable to define 
nominal SSME operating characteristics and for the 
diagnostic evaluation of system malfunction or failure. 

Extensive dynamic pressure measurements have been acquired 
at critical engine locations in the course of the SSME 
development and certification test program. Under NASA 
contract NAS 8-34343 Wyle Laboratories has developed a 
computerized data base and analysis program to characterize 
the SSME dynamic pressure environment under a wide range of 
engine operating conditions. This report summarizes the 
results of this investigation. (An interim review of the 
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program is given in Appendix A.) The objective of this 
study was threefold: 


1 . Generate a data base of dynamic pressure spectra 
obtained from SSME hot firing tests measurements. 

2. Develop data base management software to sort, 
recall and perform statistical analysis of selected 
measurements, to characterize the dynamic pressure 
environment under specified SSME operating conditions. 

3 . Develop and apply analytical techniques to 
normalize the data base and predict spectral trends 
with SSME operating condition. 


The approach employed to accomplish these objectives is 
discussed in the following sections. 

The SSME dynamic pressure data base is discussed in Section 
II. An overview of the SSME development and acceptance test 
program is given. Dynamic pressure measurements acquired 
during hot firing tests were reduced by NASA/MSFC. A brief 
description of data reduction procedures employed by MSFC is 
given and the data format provided Wyle for data base 
generation is described. Data base measurements and 
associated SSME locations are defined. Tables are included 
which summarize the SSME tests/measurements entered into the 
dynamic pressure data base. 


Section III describes the SSME dynamic pressure data base 
management system and associated analysis software. The 
original software was developed for implementation on a 
Perkin Elmer Interdata 8/32 computer system or similar 
equipment. Recently the program has been rewritten for 
application on more readily available PC/compatible systems. 
Data entry/edit routines are described along with 
subroutines for estimating statistical means, variances, 
ranges, etc. of selected measurements under user defined 
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engine operating conditions. A program listing of the data 
base management/ analysis software is included. 

A number of the dynamic pressure measurements were available 
for only a limited range of engine operating power levels. 
However, the great majority included operation at SSME rated 
power level (RPL) . In addition, an extensive data base of 
related vibration measurements at RPL has been developed by 
MSFC. Therefore ensemble average spectra were computed for 
each available measurement at RPL. These results, and their 
method of extraction are presented in Section IV. These 
summary spectra should provide a valuable reference of SSME 
dynamic pressure environments under nominal engine operating 
conditions. 

Scaling and interpolation of the dynamic pressure spectra is 
discussed in Section V. Early in the program, it was 
anticipated that dynamic pressure scaling would be performed 
through dimensional analysis based on pertinent 
flow/ geometric parameters. Several simple models were 
investigated to represent phenomena such as turbulent duct 
flow and blade wake excitation sources. A detailed review 
of data base spectra indicated, however, that the pressure 
environment is too complex to be adequately represented by a 
superposition of simple sources. Dimensional arguments were 
therefore abandoned in favor of a more generalized 
interpolation technique based on the Singular Value 
Decomposition (SVD) of a spectral matrix representing each 
measurement over the range of observed test power levels. A 
description of the SVD technique is given in Section 5-2 , 
along with an example of application to a typical data base 
measurement. The SVD technique was applied to each 
measurement representing an adequate sample size and 
available range of test power levels. In this manner a 
library of interpolated spectra was generated representing 
each measurement in 2.5% power level increments. These 
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isoplots illustrate the spectral trend of each measurement 
with varying engine power level. Finally, a software routine 
was developed to permit prediction of the pressure spectral 
density for a desired measurement at a selected engine power 
level. 

Plots representing the complete SSME dynamic pressure data 
base are included as Appendices to this report under 
separate cover. 


SECTION II 

THE SSME DYNAMIC PRESSURE DATA BASE 


2.1 Background 

The Space Shuttle Main Engines (SSMEs) are extremely 
sophisticated machines , designed to provide maximum 
performance within stringent constraints on size, weight, 
and efficiency. As a result, the SSMEs are required to 
operate under extreme temperatures with high fluid pressures 
and rotational pump speeds. Additionally, development work 
is presently in progress to uprate SSME performance. Since 
the Space Shuttle Vehicle (SSV) is man-rated, and a key to 
the cost effectiveness of the SSV concept is hardware 
reusability, system reliability is of paramount importance. 

The Orbiter vehicle main propulsion system consists of three 
SSMEs. The SSMEs are reusable, high-performance, liquid- 
propellant rocket engines with variable thrust. They are 
ignited on the ground at launch and operate in parallel, 
with approximately 500 seconds total firing duration. Each 
of the rocket engines operates at a mixture ratio (liquid 
oxygen/liquid hydrogen) of 6:1 and a chamber pressure of 
approximately 3000 psia to produce a sea-level thrust of 
375,000 pounds and a vacuum thrust of 470,000 pounds. The 
engines are presently throttleable over a thrust range of 60 
to 109 percent of the design thrust level. This provides a 
higher thrust level during liftoff and the initial ascent 
phase, and allows Orbiter acceleration to be limited to 3 
g's during the final ascent phase. The engines are gimbaled 
(+10.5 degrees for pitch and +8.5 degrees yaw) to provide 
pitch, yaw, and roll control during the Orbiter boost phase. 
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Significant to meeting performance requirements is the use 
of the staged combustion power cycle coupled with high 
combustion chamber pressures. In the SSME-staged combustion 
cycle, the propellants are partially burned at high pressure 
and relatively low temperature in the preburners, then 
completely combusted at high temperature and pressure in the 
main chamber before expanding through the high-area-ratio 
nozzle. Hydrogen fuel is used to cool all combustion 
devices in contact with high-temperature combustion 
products. An electronic engine controller automatically 
performs checkout, start, mainstage, and engine shutdown 
functions. Major components of the SSME are illustrated in 
Figure 2-1. A more detailed view of the SSME power head is 
shown in Figure 2-2. This figure provides an indication of 
the complexity of the SSME turbomachinery. 
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FIGURE 2-1. SPACE SHUTTLE MAIN ENGINE 
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FIGURE 2-2. Space Shuttle Main Engine Power Head 
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2.2 SSME Development and Acceptance Testing 


To validate system performance and ensure equipment 
reliability, the SSME and components have been and are 
presently undergoing extensive development and qualification 
tests. Testing of the engine and components is conducted at 
several NASA and contractor locations. Full scale engine 
test firings for development and flight acceptance are 
performed on two single-engine test stands at the National 
Space Technology Laboratories (NSTL) , Bay St. Louis, 
Mississippi, and at one stand operated by RocJcwell 
International near Santa Susana, California, with plans to 
include a development test stand at MSFC. In addition, main 
propulsion testing (MPT) is performed at NSTL on a stand 
designed to accommodate the Shuttle main propulsion system 
elements — the three-engine cluster, the ET, and the Orb iter 
systems . 

Testing is being performed on a continuing basis. The 
length of a given test is dependent on specific test 
objectives and may run from several seconds to over 800 
seconds. Tests are generally designed to satisfy multiple 
specific objectives, which fall into two broad categories; 
acceptance/certification firing of flight hardware and 
development testing directed towards design verification, 
performance and reliability improvement. Test operations 
are controlled by a computer called the Command and Data 
Simulator (CADS) which communicates with the engine, 
displays vital measurements for on-line observation/control 
and initiates pre- and post-test procedures. 

Approximately 250 measurements are recorded on a given test 
including wide band vibration, dynamic pressure and strain 
at critical engine locations. Some of these measurements 
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are utilized on-line as emergency cut-off indicators and all 
are recorded on magnetic tape for subsequent analysis and 
evaluation. Limited SSME measurements are recorded on 
magnetic tape during SSV flights for evaluation with orbiter 
return. Typical dynamic measurements obtained during SSME 
operation are illustrated in Figure 2-3. 

2.3 Data Base Measurements. Format and Locations 

The dynamic pressure data provided by MSFC was received in 
three formats for defining each measurement: 

1. Engine Power Level as a function of time. 

2 . Root-Mean-Square pressure as a function of time 
(0.4 second averaging time) 

3. Mean-square density spectra (PSD's) defining the 
frequency decomposition of each measurement at 
selected test conditions. 

The spectral data were defined over 1-5 kHz, 1-20 kHz, and 
1-40 kHz frequency bands. The great majority of data were 
reduced over a 1-20 kHz frequency range. The spectral data 
were reduced on a B&K wave analyzer yielding 400 spectral 
estimates over the range of analysis. This proved to be a 
consideration in the data base entry from continuous analog 
plots provided. 

The dynamic pressure measurements are defined in terms of a 
number of acronyms for efficiency in computerized log and 
test description. A glossary of these abbreviations is 
given in Table 2-1. Table 2-2 defines the measurements and 
engine locations utilized in data base development. The 
SSME test measurement locations are illustrated in Figure 2- 
3 . A summary of test measurements thus far entered into the 
data base is given in Tables 2-3, 2-4, and 2-5. 
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ENGINE FUEL IN PR, TEW 
( 0F3F9, 0F3F8D ) 


INSTRUMENTATION 


OXIDIZER 


LPOP RAD, AX ACC 
(AZ2AA, ZOA, 2CA, 1AA ) 



view A-A 


NOZZLE AFT NAN PR 1 TEMP 
( NOZ 2 ) 




VIEW C-C 
tc. INJECTOR 


□ OPB PC 
(OGIC) 


viewB-B 



,0X1 D PB PC 
(OGIB) 


SFiN 

o 



\ 

OPB PC 
(OGIA) 


O 

□ 


□ 



G 



O 


Hi irwi imcih vtatminM 
nst larngMatitiOM n Misses 
nun iittinuiuiM 


FLIGHT INSTRUMENTATION (CADS) 
APPENDED DATA WORDS (CADS) 
BECKMAN DATA (BKM) 

BECKMAN NFD DATA 

HIGH FREQUENCY ANALOG DATA 


2/15/79 


FIGURE 2-3. 


TYPICAL SSME MEASUREMENT LAYOUT 
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TABLE 2-1. GLOSSARY OF ABBREVIATIONS USED 
IN MEASUREMENT DESCRIPTION 


LPFP 

Low Pressure Fuel Pump 

IN 

Inlet 

PR 

Pressure 

LPOP 

Low Pressure Oxider Pump 

DS 

Discharge 

HPOP 

High Pressure Oxidizer Pump 

BAL 

Balance 

CAV 

Cavity 

LOX 

Liquid Oxygen 

DISC 

Discharge 

MCC 

Main Combustion Chamber 

PBP 

Prebumer Pump (Oxidizer) 

SUP 

Supply 

FPB 

Fuel Prebumer 

LPOT 

Low Pressure Oxidizer Turbine 

TURB 

Turbine 

DR 

Drive 

HX 

High Pressure 

OPB 

Oxidizer Prebumer 

HEX 

High Pressure Lox 

HPFP 

High Pressure Fuel Pump 

MAN 

Manifold 

INJ 

In j ector 

PC 

Chamber Pressure 

DC 

(Direct Current) Indicates that 
measurement includes mean pressure 

AUX 

Auxiliary 

LX 

Lox 

HI 

High Pressure 

OX 

Oxidizer 
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TABLE 2-2. DYNAMIC PRESSURE MEASUREMENT 
DESCRIPTIONS AND SSME LOCATION (FROM BECKMAN TEST LOG) 


MEASUREMENT DESCRIPTION 

LPFP in PR 

LPOP DS PR 

HPOP BAL CAV PR1 

HPOP DS PR 

LOX IN DUCT PR2 

HPOP DISC PR 

MCC HOT GAS IN PR 

PBP DS PR 

CHAMBER PR 

LPFP DS PR 

PBP FUEL SUP PR 

FPB FUEL MAN PR 

LPOT TURB DR PR 

LOX HX OUT PR 

OPB PC PR 

HEX OUTLET 

HPFP DS PR 

HOT GAS MAN PR3 

HPFP BAL CAV PR 

MCC FUEL INJ PR 

FPB PC 

MCC FUEL INJ PR DC 
OPB PC 
HPOP IN PR 
AUX LX I PR 
HI LX I PR 
HOT GAS MAN PR2 
MCC OX INJ PR 
OPB PC PR DC 
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The management, analysis and interpolation of these data to 
characterize the SSME dynamic pressure environment is 
described in the following sections. 
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TABLE 2-3. DATA PLOT LISTING FOR TEST STAND A1 


Test 

290 


T 

293 


294 


Channel Descriptor 
(Time History or PSD) 


Maximum Time 

Frequency Slice 

(Hz) (Sec) 


MMC HOT GAS IN PR 


HPOP DS PR 

I 

PBP DS PR 

T 

OPB PC PR DC 

v 

HPOP DS PR 


PBP DS PR 


T 


OPB 


PL 


40.000 25 
55 
75 

195 

315 

,, 435 

500 

20.000 10 
200 
400 
600 
800 

10 

200 

400 

600 

800 

10 

200 

400 

, r 600 

" 800 

20,000 10 

40 
100 
300 
500 
645 
660 
10 
40 
100 
300 
500 
645 
660 
10 
40 
100 
300 
500 

,, 645 

* 660 
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DATA PLOT LISTING FOR TEST STAND A1 (Continued) 


Test 

294 


295 


297 

I 

299 


300 

I 


Channel Descriptor 
(Time History or PSD) 

HPOP IN PR 


HPOP DS PR 


» 

PBP DS PR 


HPOP IN PR 


f 

OPB PC PR DC 


HPOP DS PR 
PBP DS PR 
OPB PC PR DC 
HPOP DS PR 
PBP DS PR 
HPOP IN PR 
AUX LX I PR 
OPB PC PR DC 
PBP DS PR 
HPOP IN PR 
OPB PC PR DC 


Maximum Time 

Frequency Slice 

(Hz) (Sec) 


20,000 10 


40 

100 

300 

500 

645 

660 

20 

45 

100 

300 

450 

500 

20 

45 

100 

300 

450 

500 

20 

45 

100 

300 

450 

500 

20 

45 

100 

300 

450 

500 

6 

I 

6 
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DATA PLOT LISTING FOR TEST STAND A1 (Continued) 


Test 

301 


303 


304 


305 

I 


Channel Descriptor 
(Time History or PSD) 

HPOP DS PR 

I 

PBP DS ER 

\ 

OPB PC PR DC 

I • 

HPOP DS PR 

I 

PBP DS P fc 

I 

HPOP IN PR 

\ 

OPB PC PR DC 

I 

? 

HPOP DS PR 

I 

PBP DS PR 

I 

LPFP DS PR 

I 

FPB FUEL MAN PR 

I 

HPOP IN PR 

I 

HI LX I PR 

I 

HPOP DS PR 


Maximum Time 

Frequency Slice 

(Hz) (Sec) 


20,000 100 


250 

400 

100 

250 

400 

100 

250 

400 

5 

30 

45 

5 

30 

45 

5 

30 

45 

5 

30 

45 

5 

50 

75 

5 

50 

75 

5 

50 

75 

5 

50 

75 

5 

50 

75 

5 

50 

75 

5 

25 

50 

75 
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Test 

305 


306 


307 

I 


DATA PLOT LISTING FOR TEST STAND A1 (Continued) 


Channel Descriptor 
(Time History or PSD) 


Maximum Time 

Frequency Slice 

(Hz) (Sec) 


PBP DS PR 


FPB FUEL MAN PR 


HPOP IN PR 


HI LX I PR 


HPOP DS PR 


MCC HOT GAS IN PR 


PBP DS PR 


FPB FUEL MAN PR 


HPOP IN PR 


HI LX I PR 


HPOP DS PR 


20,000 


40,000 


20,000 


5 

25 

50 

75 

5 

25 

50 

75 

5 

25 

50 

75 

5 

25 

50 

75 

5 

6 
35 
70 

5 

20 

50 

75 

5 

35 

70 

5 

20 

50 

75 

5 

20 

50 

75 

5 

20 

50 

75 

5 

20 

45 

70 
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DATA PLOT LISTING FOR TEST STAND A1 (Continued) 




Maximum 

Time 


Channel Descriptor 

Frequency 

Slice 

Test 

(Time History or PSD) 

(Hz) 

(Sec) 

307 

MCC HOT GAS IN PR 

40, 

000 

5 



1 

1 


20 



1 

1 


45 



f 

PBP DS PR 

1 

f 

70 



20, 

000 

5 



1 



20 



1 



45 



t 



70 



FPB FUEL MAN PR 



5 



1 



20 



] 



45 



t 



70 



HPOP IN PR 



5 



1 



20 



1 



45 



f 



70 



HI LX I PR 



5 



i 



20 


f 

T 



70 

309 

HPOP. DS PR 



6 



t 



30 



MCC HOT. GAS IN PR 



6 



t 



24 



PBP DS PR 



6 



f 



30 



HPOP IN PR 



6 



1 



24 

1 


HI LX I PR 

- 


6 

1 


f 



24 

310 

HPOP DS PR 



6 

1 


t 



30 



MCC HOT GAS IN PR 



6 

i 


f 



30 



PBP DS PR 



6 



f 


1 

30 



HPOP IN PR 

2, 

ooo 

6 



♦ 

1 

f 

30 



HI LX I PR 

20, 

ooo 

6 



f 



30 

3: 

.1 

HPOP.DS PR 



6 



♦ 



30 



MCC HOT GAS IN PR 



6 

i 


1 



30 

1 


PBP PS PR 



6 

1 


t 



30 


2-15 


DATA PLOT LISTING FOR TEST STAND A1 (Continued) 



31 


316 


317 


Channel Descriptor 
(Time History or PSD) 

HPOP IN PR 

t 

HI LX I PR 

▼ 

mcc hot Gas in pr 

T 

HPOP IN PR 

1 

HI LX I PR 

MCC HOT GAS IN PR 

I 

HPOP IN PR 

T * 

HI LX I PR 

T 

HPOP DS PR 

▼ 

MCC HOT GAS IN PR 


t 

PBP DS PR 

T 

HPOP IN PR 

1 

HI LX I PR 

▼ 

HPOP DS PR 

f 

MCC HOT GAS IN PR 

▼ 

PBP DS PR 

T 

HPOP DS PR 

f 

MCC HOT GAS IN PR 


PBP DS PR 

f 

OPB PC PR DC 

f 

HPOP DS PR 


Maximum 

Frequency 

(Hz) 


Time 

Slice 

(Sec) 


20,000 


40.000 

I 

20.000 


40.000 

T 

20.000 


t 


6 

30 

6 

30 

6 

40 

6 

40 

6 

4 

30 

4 

30 

3 
30 

6 

30 

6 

30 

6 

30 

6 

30 

6 

30 

6 

30 

4 
30 
30 
60 

4 

30 

4 

11 

4 

11 

4 

11 

4 

11 

4 

11 

4 

11 
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DATA PLOT LISTING FOR TEST STAND A1 (Continued) 


Maximum Time 

Channel Descriptor Frequency Slice 

Test (Time History or PSD) (Hz) (Sec) 


317 


T 

319 


MCC HOT GAS IN PR 


t 

PBP DS -PR 

? 

OPB PC PR DC 

t 

HPOP DS PR 


321 


MCC HOT GAS IN PR 

I • 

PBP DS PR 

I 

MCC IN PR 

I 

OPB PG PR DC 
MCC HOT GAS IN PR 


FPB PC 


OPB PC PR DC 


20,000 

▼ 

40.000 

T 

20.000 


40.000 

I 

20.000 


40,000 


t 

20,000 


4 

11 

4 

11 

4 

•11 

4 

11 

6 

30 

92 

6 

30 

92 

6 

30 

92 

12 

60 

180 

6 

30 

92 

6 

30 

90 

100 

270 

6 

30 

90 

100 

270 

6 

30 

90 

100 

270 

6 

30 

90 

100 

270 
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DATA PLOT LISTING FOR TEST STAND A1 (Continued) 


Maximum Time 

Channel Descriptor Frequency Slice 

Test (Time History or PSD) (Hz) (Sec) 


322 HPOP DS PR 


1 

MCC HOT GAS IN PR 


I 

PBP DS PR 


FPB PC 


! 

OPB PC PR DC 


t 


20,000 6 
130 
170 
190 
210 
230 
260 
280 
6 

130 

150 

170 

190 

210 

230 

260 

280 

6 

130 

150 

170 

190 

210 

230 

264 

284 

6 

130 

150 

170 

190 

210 

230 

260 

280 

6 

130 

150 

170 

190 

210 

230 

,, 260 

' 280 
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DATA PLOT LISTING FOR TEST STAND A1 (Continued) 


Maximum Time 

Channel Descriptor Frequency Slice 

Test (Time History or PSD) (Hz) (Sec) 


384 

385 

388 

389 

390 

391 

393 

394 

395 

397 

398 

400 

401 

402 

403 

404 


THRUST 


405 

406 

407 


HPOP BAL CAV PR2 

I 

HPOP BAL CAV PR1 

I 

OPB PC PR 

I 

HPFP BAL CAV PR 

I 

FPB PC 

I 

THRUST 

I 

HPOP BAL CAV PR2 


HPOP BAL CAV PR1 


20,000 20 
100 
220 
20 
100 
220 
20 
100 
220 
20 
100 
220 
20 

If 100 

" 220 


20,000 20 

45 
75 
200 
400 
430 
20 
45 
75 
200 

,i 400 

430 
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DATA PLOT LISTING FOR TEST STAND A1 (Continued) 


Maximum Time 

Channel Descriptor Frequency Slice 

Test (Time History or PSD) (Hz) (Sec) 


407 


FPB PC 


f 

408 

409 


OPB PC PR 


T 

HPFP BAL CAV PR 


THRUST 

1 

HPOP BAL CAV PR2 


HPOP BAL CAV PR1 


T 

OPB PC PR 


HPFP BAL CAV PR 


! 


20,000 20 
45 
75 
200 
400 
430 
20 
45 
75 
200 
400 
430 
20 
45 
75 
200 
400 
430 


20,000 15 

240 
500 
525 
705 
725 
15 
240 
500 
525 
705 
725 
15 
240 
500 
525 
705 
725 
15 
240 
500 
525 

,, 705 

725 
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Test 

409 

410 

413 

414 


DATA PLOT LISTING FOR TEST STAND A1 (Concluded) 


Maximum Time 

Channel Descriptor Frequency Slice 

(Time History or PSD) (Hz) (Sec) 

FPB PC 20,000 7 

240 

500 

525 

m ' ,, 705 

' T 725 

THRUST 

I . 
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TABLE 2-4. DATA PLOT LISTING FOR TEST STAND A2 


Channel Descriptor 
Test (Time History or PSD) 

214 MMC FUEL INJ PR 


215 

t I 

251 THRUST 

252 

253 

254 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

276 

277 
279 

283 

284 
303 

306 

307 

310 

311 

313 „ 

314 » 


Maximum Time 

Frequency Slice 

(Hz) (Sec) 


40,000 25 


50 

70 

150 

300 

450 

500 

30 

50 

70 
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TABLE 2-5. DATA PLOT LISTING FOR TEST STAND A3 


Maximum Time 

Channel Descriptor Frequency Slice 

Test (Time History or PSD) (Hz) (Sec) 


) 


\ 


t 

( 


047 THRUST 

048 

049 

050 

052 

053 
066 

067 

068 

069 

070 

071 

072 

073 

074 

075 

076 

077 

078 

079 

087 

088 

089 

090 

091 

092 

093 

094 

095 

096 

097 

098 

099 ,, 

100 

LPFP IN PR 
LPOP DS PR 
HPOP BAL CAV PR2 
HPOP BAL CAV PR1 
LPFP IN PR 


5,000 

\ 


10 

40 

85 

95 


i 
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DATA PLOT LISTING FOR TEST STAND A3 (Continued) 


Maximum Time 

Channel Descriptor Frequency Slice 

Test (Time History or PSD) (Hz) (Sec) 


100 LPFP IN PR 


t 

LPOP DS PR 


1 

HPOP BAL CAV PR2 


V* 

HPOP BAL CAV PR1 


101 THRUST 

LPFP IN PR 
LPOP DS PR 
HPOP BAL CAV PR2 
LPFP IN PR 


10 


LPOP DS PR 


HPOP BAL CAV PR2 


THRUST 
LPFP IN PR 
HPOP BAL CAV PR2 
HPOP DS PR 
LPFP IN PR 

I 


20,000 10 
40 
85 
95 
10 
40 
85 
95 
10 
40 
85 
95 
10 
40 

,, 85 

f 95 


5,000 5 

I 30 

1 80 

f 95 

20,000 5 

30 
80 
95 
5 
30 
80 
95 
5 
30 

ii so 

95 


20,000 15 

35 
60 
80 
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DATA PLOT LISTING FOR TEST STAND A3 (Continued) 


Test 

102 


103 


Maximum Time 

Channel Descriptor Frequency Slice 

(Time History or PSD) (Hz) (Sec) 


LPFP IN PR 

• I 

HPOP BAL CAV PR2 


HPOP DS PR 


LOX IN DUCT PR2 


THRUST 
LPFP IN PR 
LPOP DS PR 
HPOP BAL CAV PR2 
HPOP BAL CAV PR1 
HPOP DISC PR 
LPFP IN PR 


20,000 


5,000 

20,000 


100 

200 

280 

15 

35 

60 

80 

100 

150 

200 

280 

15 

35 

60 

80 

100 

200 

280 

15 

35 

60 

200 

280 

15 

35 

60 

80 

100 

200 

280 


5,000 


f 


20,000 

* 


5 

15 

20 

26 

30 

5 

15 

20 
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DATA PLOT LISTING FOR TEST STAND A3 (Continued) 


Test 

103 


104 


105 

T 


Channel Descriptor 
(Time History or PSD) 

LPFP IN PR 

▼ 

LPOP DS PR 


HPOP BAL CAV PR2 


HPOP DISC PR 


f 

MCC HOT GAS IN PR 


MMC LONG ACCEL 1 
THRUST 
LPFP IN PR 
LPOP DS PR 
HPOP BAL CAV PR2 
HPOP BAL CAV PR1 
HPOP DISC PR 
MMC HOT GAS IN PR 
LPFP IN PR 
LPOP DS PR 
HPOP BAL CAV PR2 
HPOP DISC PR 
MMC HOT GAS IN PR 
? 

PBP DS PR 
THRUST 
LPFP IN PR 


Maximum Time 

Frequency Slice 

(Hz) (Sec) 


20,000 26 


30 

5 

15 

20 

26 

30 

5 

15 

20 

26 

30 

5 

15 

20 

26 

30 

5 

15 

20 

26 

30 

40,000 5 

15 
20 
26 
30 
5 


20,000 6 


t 

40.000 

20.000 


2-26 


DATA PLOT LISTING FOR TEST STAND A3 (Continued) 


Maximum Time 

Channel Descriptor Frequency Slice 

Test (Time History or PSD) (Hz) (Sec) 


105 LPOP DS PR 
HPOP BAL CAV PR2 
HPOP BAL CAV PR1 

HPOP DS PR 
* HPOP DISC PR 
MMC HOT GAS IN PR 
PBP DS PR 
LPFP IN PR 
LPOP DS PR 
HPOP BAL CAV PR2 
HPOP BAL CAV PR1 
HPOP DS PR 
HPOP DISC PR 
MMC HOT GAS IN PR 
▼ 

, , PBP DS PR 

CHAMBER PR 

106 THRUST 
LPFP IN PR 
LPOP DS PR 

HPOP BAL CAV PR2 
HPOP BAL CAV PR1 
HPOP DS PR 
HPOP DISC PR 
MCC HOT GAS IN PR 
PBP DS PR 
LPFP IN PR 


t 

LPOP DS PR 

▼ 

HPOP BAL CAV PR2 

? 

HPOP BAL CAV PR1 
HPOP DS PR 

▼ 

HPOP DISC PR 

T 

, , PBP DS PR 

T I 

107 THRUST 


20,000 6 


40.000 

20.000 

? 


5,000 6 

T 11 

20,000 6 

11 
6 

11 

6 

11 

11 

6 

11 

6 

11 

< f 6 

11 
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DATA PLOT LISTING FOR TEST STAND A3 (Continued) 


Test 

108 


109 


T 

110 



Maximum 

Time 

Channel Descriptor 

Frequency 

Slice 

tTime History or PSD) 

(Hz) 

(Sec) 

THRUST 



LPFP IN PR 

5,000 

6 

1 

f 

12 

1 

20,000 

6 

V 


12 

LPOP DS PR 


6 

» 


12 

HPOP DS PR 


6 

? 


12 

HPOP DISC PR 


6 

? 


12 

MCC HOT GAS IN PR 


6 

1 * 

40,000 

6 

t 

I 

12 

PBP DS PR 

20,000 

6 

1 

T 

12 

THRUST 




LPFP IN PR 
LPOP DS PR 
HPOP BAL CAV PR2 
HPOP BAL CAV PR1 
HPOP DS PR 
HPOP DISC PR 
MCC HOT GAS IN PR 
PBP DS PR 
LPFP IN PR 
I 


20,000 


6 

14 


LPOP DS PR 

T 

THRUST 
LPFP DS PR 
LPOP DS PR 
HPOP BAL CAV PR2 
HPOP BAL CAV.PR1 
CHAMBER PR 
PBP FUEL SUP PR 
FPB FUEL MAN PR 
LPFP DS PR 


f 

HPOP BAL CAV PR2 

t 


I 


6 

14 


20,000 10 

30 
150 
300 
400 
10 

,i 30 

f 150 
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DATA PLOT LISTING FOR TEST STAND A3 (Continued) 


Maximum Time 
Channel Descriptor Frequency Slice 

Test (Time History or PSD) (Hz) (Sec) 


110 HPOP BAL CAV PR2 

T 

HPOP DS PR 


MCC HOT GAS IN PR 


I ' 

111 THRUST 

LPFP DS PR 
LPOP DS PR 
HPOP BAL CAV PR2 
HPOP BAL CAV PR1 
HPOP DS PR 
BPB DS PR 
CHAMBER PR 
PBP FUEL SUP PR 
FPB FUEL MAN PR 
LOX HX OUT PR 
LPFP DS PR 


LPOP DS PR 


20,000 


T 


40,000 

I 


300 

400 

10 

30 

150 

300 

10 

30 

150 

10 

30 

150 


20,000 15 

40 
50 
60 
70 
80 
90 
100 
120 
200 
280 
15 
40 
50 
60 
70 
80 
90 
100 
120 

,, 200 

280 
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DATA PLOT LISTING FOR TEST STAND A3 (Continued) 


Test 

111 


Channel Descriptor 
(Time History or PSD) 

HPOP BAL CAV PR2 


HPOP BAL CAV PR1 


f 

HPOP DS PR 


PBP DS PR 


Maximum Time 

Frequency Slice 

(Hz) (Sec) 


20,000 15 


40 

50 

60 

70 

80 

90 

100 

120 

200 

280 

15 

40 

50 

60 

70 

80 

90 

100 

120 

200 

280 

15 

40 

50 

60 

70 

80 

90 

100 

120 

200 

280 

15 

40 

50 

60 

70 

80 

90 

100 

120 

200 

280 
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DATA PLOT LISTING FOR TEST STAND A3 (Continued) 


Test 

111 


Channel Descriptor 
(Time History or PSD) 


Maximum Time 

Frequency Slice 

(Hz) (Sec) 


PBP FUEL SUP PR 


T 

LPOT TURB DR PR 


T 

LOX HX OUT PR 


20,000 


5,000 


20,000 


5,000 


20,000 


15 

40 

50 

60 

70 

80 

90 

100 

120 

200 

280 

15 

40 

80 

90 

100 

120 

200 

1-5 

40 

50 

60 

70 

80 

90 

100 

120 

200 

280 

15 

40 

50 

60 

70 

80 

90 

100 

120 

200 

280 

40 

50 

60 

70 
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DATA PLOT LISTING FOR TEST STAND A3 (Continued) 


Test 

Channel Descriptor 
(Time History or PSD) 

Maximum 

Frequency 

(Hz) 

Time 

Slice 

(Sec) 

111 

LOX HX OUT PR 

20,000 

80 
90 
100 
' 120 
200 
280 
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DATA PLOT LISTING FOE TEST STAND A3 (Continued) 


Maximum Time 

Channel Descriptor Frequency Slice 

Test (Time History or PSD) (Hz) (Sec) 


112 THRUST 

LPOP DS PR 
HPOP BAL CAV PR2 
HPOP BAL CAV PR1 
HPOT TURB DR PR 
LPFP DS PR 

* 

LPOP DS PR 

i 

HPOP BAL CAV PR2 

i ■ 

MCC HOT GAS IN PR 


» 

LPOT TURB DR PR 


T 

OPB PC PR 


HEX OUTLET 

* 


20.000 30 
150 
280 

30 

150 

280 

30 

150 

280 

30 

, 150 

T 280 

40.000 30 

1 150 

T 280 

5,000 30 

I 150 

f 280 

20.000 30 
150 
280 

30 

150 

280 

30 

w 150 

T 280 


1 


3 THRUST 

LPOP DS PR 
HPOP BAL CAV PR2 
HPOP BAL CAV PR1 
HPOP DS PR 
MCC HOT GAS IN PR 
PBP DS PR 


LPOP DS PR 

20,000 

30 

1 

i 

150 

t 

t 

280 

, HPOP BAL CAV PR2 

5,000 

150 

f 

T 

280 


2-33 


DATA PLOT LISTING FOR TEST STAND A3 (Continued) 


Test 

113 


114 

* 


Maximum Time 

Channel Descriptor Frequency Slice 

(Time History or PSD) (Hz) (Sec) 


HPOP BAL CAV PR2 


HPOP DS PR 


MCC HOT GAS IN PR 


f 

PBP DS PR 


LPFP DS PR 

I 

LPOT TURB DR PR 


OPB PC PR 


MCC LONG ACCEL 1 
OPB PC PR DC 


THRUST 
LPOP DS PR 
HPOP BAL CAV PR2 


20,000 

I 

5,000 

I 

20,000 


40.000 
5,000 

I 

20.000 


5,000 

I 

20,000 

I 

5,000 

I 

20,000 

40.000 
5,000 

{ 

20.000 

* 


30 

150 

280 

30 

150 

280 

30 

150 

280 

30 

150 

280 

30 

150 

30 

150 

280 

30 

150 

280 

30 

150 

280 

30 

150 

280 

30 

150 

280 

30 

150 

280 

30 

150 

280 

280 

30 

150 

280 

30 

150 

280 
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DATA PLOT LISTING FOR TEST STAND A3 (Continued)i 


Test 

114 


115 


Channel Descriptor 
(Time History or PSD) 


HPOP BAL CAV PR1 
HPOP DS PR 
MCC HOT GAS IN PR 
PBP DS PR 
LPOT TURB DR PR 
LPOP DS PR 

I 

HPOP BAL CAV PR2 

I 

HPOP DS PR 

i ' 

MCC HOT GAS IN PR 

I 

PBP DS PR 

l 

LPFP DS PR 

I 

LPOT TURB DR PR 


f 

OPB PC PR 

I 

HPOP DS PR 

* 

THRUST 
L’POP DS PR 
HPOP BAL CAV PR2 
HPOP BAL CAV PR1 
MCC HOT GAS IN PR 
PBP DS PR 
HPFP DS PR 
LPOP DS PR 

* 


Maximum Time 

Frequency Slice 

(Hz) (Sec) 


20,000 30 

150 
280 
30 
150 
280 
30 
150 
280 
30 
150 
280 
30 
150 
280 
30 
150 
280 

5,000 30 

I 150 

T 280 

20,000 30 

150 
280 
30 
150 
280 
30 

,, 150 

' 280 


20,000 30 

I 150 

f 280 
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Test 

115 


116 


DATA PLOT LISTING FOR TEST STAND A3 (Continued) 


Maximum Time 
Channel Descriptor Frequency Slice 

(Time History or PSD) (Hz) (Sec) 


HPOP BAL CAV PR2 


MCC HOT GAS IN PR 


PBP D.S PR 


LPFP DS PR 

I 

LPOT TURB DR PR 


OPB PC PR 

1 

HPFP DS PR 


I 

HOT GAS MAN PR3 

t 

OPB PC PR DC 

I 


THRUST 
LPOP DS PR 
HPOP BAL CAV PR2 
HPOP BAL CAV PR1 
HPOP DS PR 
MCC HOT GAS IN PR 
PBP DS PR 
OPB PC PR 
HPFP DS PR 
HOT GAS MAN PR3 
HPFP BAL CAV PR 


5,000 30 

1 150 

f 280 

20,000 30 

150 
280 
30 

w 150 

' 280 

40.000 30 

I 150 

f 280 

20.000 30 
150 
280 

30 

150 

280 

5,000 30 

I 150 

f • 280 
20,000 30 

150 
280 
30 
150 
30 
150 
280 
30 
150 
280 
30 

u 150 

T 280 
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Test 

116 


DATA PLOT LISTING FOR TEST STAND A3 (Continued) 


Channel Descriptor 
(Time History or PSD) 


LPOP DS PR 


HPOP BAL CAV PR2 


HPOP DS> PR 

I 

MCC HOT GAS IN PR 


Maximum 

Time 

Frequency 

Slice 

(Hz) 

(Sec) 

5,000 

30 

1 

150 

T 

280 

20,000 

30 

[ 

150 

t 

280 

5,000 

30 

1 

150 

♦ 

280 

20,000 

30 


PBP DS PR 


40.000 

5,000 

I 

20.000 


LPFP DS PR 


LPOT TURB DR PR 


OPB PC PR 
HPFP DS PR 

I 

HOT GAS MAN PR3 


5,000 

I 

20,000 
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DATA PLOT LISTING FOR TEST STAND A3 (Continued) 


Maximum Time 

Channel Descriptor Frequency Slice 

Test (Time History or PSD) (Hz) (Sec) 


116 . HPFP BAL CAV PR 

I 

OPB PC PR DC 

I 1 ' 

117 THRUST 
LPOP DS PR 
HPOP DS PR 

MCC HOT GAS IN PR 
PBP DS PR 
LPFP Dg PR 
LPOT TURB DR PR 
OPB PC PR 
HPFP DS PR 
HPFP BAL CAV PR 
LPOP DS PR 


HPOP DS PR 

I 

MCC HOT GAS IN PR 


V 

PBP DS PR 

* 

LPFP DS PR 

I 

LPOT TURB DR PR 

I 

OPB PC PR 

I 


20,000 30 

150 
280 
4 
30 

' 150 

280 


5,000 30 

I 150 

f 297 

20.000 30 
150 
297 

30 

150 

297 

30 

if 1^0 

" 297 

40.000 30 

I 150 

f 280 

20.000 30 
150 
297 

30 

150 

297 

30 

150 

297 

30 

ii , 150 

' 297 
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DATA PLOT LISTING FOR TEST STAND A3 (Continued) 


Test 

117 


T 

118 


Channel Descriptor 
(Time History or PSD) 

HPFP DS PR 

I 

HPFP BAL CAV PR 

- I ’ 

OPB PC PR DC 

* . 

THRUST 
LPOP DS PR 
HPOP BAL CAV PR1 
PBP DS PR 
LPOP DS PR 

I 

MCC HOT GAS IN PR 


PBP DS PR 

I 

HPFP BAL CAV PR 

I 

OPB PC PR DC 

I 


Maximum Time 

Frequency Slice 

(Hz) (Sec) 


20,000 30 


150 

297 

30 

150 

297 

30 

150 

297 


20.000 30 
150 
280 

30 

i f 150 

280 

40.000 30 

I 150 

T 280 

20.000 30 
150 
280 

30 

150 

280 

30 

. r 150 

280 


2-39 


Test 

120 


t 

121 


9 


DATA PLOT LISTING FOR TEST STAND A3 (Continued) 


Maximum Time 

Channel Descriptor Frequency Slice 

(Time History or PSD) (Hz) (Sec) 


THRUST 

HPOP BAL CAV PR2 
HPOP BAL CAV PR1 
HPOP DS PR 
MCC HOT GAS IN PR 
PBP DS PR 
LPOP DS PR 

I 

HPOP DS PR 

I 

MCC HOT GAS IN PR 


PBP DS PR 

I. 

HPFP BAL CAV PR 

I 

OPB PC PR DC 

I 

THRUST 
LPOP DS PR 

T 

HPOP DS PR 

I 

MCC HOT GAS IN PR 


HOT GAS MAN PR2 

T 

HPFP BAL CAV PR 

T 

OPB PC PR DC 

▼ 


20,000 30 
100 
200 
30 
100 
200 
30 
100 

* 200 

40.000 30 

I 100 

♦ 200 

20.000 30 
100 
200 

30 

100 

200 

30 

100 

" 200 


20,000 8 

16 
8 

16 

8 

" 16 

40.000 8 

f 16 

20.000 8 

16 

8 

16 

w 8 

? 16 
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Test 

123 

124 

125 


126 

127 

I 

128 

1 

129 


130 


DATA PLOT LISTING FOR TEST STAND A3 (Continued) 


Maximum Time 

Channel Descriptor Frequency Slice 

(Time History or PSD) (Hz) (Sec) 


THRUST 

I 

MCC FUEL INJ PR 20,000 6 


p 95 

125 

155 

185 

215 

245 

I 265 

♦ 285 

„ 40,000 6 

10 

50 

95 

125 

155 

185 

215 

245 

265 

T 285 

THRUST 

▼ 


FPB PC 

20,000 

6 

1 

r 

1 

r 

11 

THRUST 




FPB PC 

20,000 

6 

OPB PC PR DC 

1 

r 

f 

THRUST 




FPB PC 

20,000 

6 

t 

! 

f 

11 

MCC FUEL INJ PR 

40,000 

5 

\ 

? 

1 

r 

11 

OPB PC PR DC 

20,000 

6 

1 

r 

1 

r 

11 

THRUST 




FPB PC 

20,000 

6 


T 



11 

MCC FUEL INJ PR 



6 





11 



40,000 

6 



f 

11 


2-41 


DATA PLOT LISTING FOR TEST STAND A3 (Continued) 


Channel Descriptor 
Test (Time History or PSD) 


Maximum Time 

Frequency Slice 

(Hz) (Sec) 


130 

▼ 

131 


OPB PC PR DC 

T 

THRUST 

MCC FUEL INJ PR 

I ' 

FPB PC 


OPB PC PR.DC 

I 


132 


THRUST 

MCC FUEL INJ PR 
OPB PC PR DC 
FPB PC 


MCC FUEL INJ PR 


OPB PC PR DC 


133 


THRUST 
FPB PC 


20,000 

1 

20,000 


40, 


000 


20,000 

♦ 


6 

11 

30 

70 

100 

30 

70 

100 

30 

70 

100 

30 

70 

100 


20,000 


V 


40,000 


20,000 


10 

90 

150 

260 

291 

10 

90 

150 

260 

291 

10 

90 

150 

260 

291 

10 

90 

150 

260 

291 


20,000 6 
71 
76 
150 

" 280 
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DATA PLOT LISTING FOR TEST STAND A3 (Continued) 


Maximum Time 

Channel Descriptor Frequency Slice 

Test (Time History or PSD) (Hz) (Sec) 


133 


OPB PC PR DC 


134 


THRUST 

MCC FUEL INJ PR 


FPB PC 


MCC FUEL INJ PR DC 


OPB PC PR DC 


135 THRUST 

I MCC FUEL INJ PR 


20,000 6 
71 
76 
150 

f 280 


20,000 


40,000 


20,000 


40,000 


20,000 


4 

30 

45 

150 

280 

4 

30 

45 

150 

200 

4 

30 

55 

150 

280 

4 

30 

45 

150 

280 

4 

30 

45 

150 

280 

4 

30 

45 

150 

280 


20,000 6 

15 
45 
75 

11 150 

40,000 6 

15 
45 
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DATA PLOT LISTING FOR TEST STAND A3 (Continued) 


Test 

135 


t 

136 


138 

139 

140 

141 

142 

143 

144 

145 

146 


Channel Descriptor 
(Time History or PSD) 

MCC FUEL INJ PR 

T 

MCC FUEL INJ PR DC 


OPB PC PR DC 


THRUST 

MCC FUEL INJ PR 


MCC FUEL INJ PR DC 


OPB PC PR DC 


THRUST 


Maximum 

Frequency 

(Hz) 

40.000 

V 

20.000 


40,000 


20,000 


20,000 


40,000 


20,000 


Time 

Slice 

(Sec) 

75 

150 

6 

15 

45 

.75 

150 

6 

15 

45 

75 

150 

6 

15 

45 

75 

150 

4 

25 

50 

150 

280 

4 

25 

50 

150 

280 

4 

25 

50 

150 

280 

25 

50 

150 

280 


2-44 


DATA PLOT LISTING FOR TEST STAND A3 (Concluded) 


Channel Descriptor 
Test (Time History or PSD) 


Maximum 

Frequency 

(Hz) 


147 

148 

150 

151 

152 

153 

154 

155 

156 

157 

158 
161 

164 

165 

166 

167 

168 

170 

171 

172 

173 

174 

175 

177 

178 

179 

180 
181 
182 

183 

184 

185 

186 

187 

188 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 


THRUST 


Time 

Slice 

(Sec) 
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Section III 


THE SSME DYNAMIC PRESSURE DATA BASE 
MANAGEMENT AND ANALYSIS SYSTEM 

A primary objective of this study was to develop a data base 
and data base management system for SSME dynamic pressure 
measurements provided by MSFC. The data was obtained in 
the form of time histories (power level and r.m.s. pressure) 
and power spectra (PSDs) . The data base management and 
analysis system was originally designed for use on an 
Interdata 8/32 computer system and other computer systems 
utilizing similar file name structures. The data 
base/management system has recently been converted for 
application on P.C. network systems. The data management 
system is designed to store test data in a manner that will 
provide rapid and easy access to the data by the display and 
analysis programs. The system can accept time history data 
and PSD or frequency domain data. Each time history and PSD 
is stored in a separate data file. The data base management 
system creates and names these files when the data is 
entered. The system also makes entries in several directory 
files so that the data can be quickly accessed, sorted, 
statistically summarized and or plotted as required by the 
user. 

Figures 3-1 through 3-3 show sample plots of data available 
in the data base. Separate volumes of plots of each data 
file in the data base have been prepared, and are provided 
in appendices to this report. 

This section describes the data base management system, the 
software that has been developed to enter data into the data 
base system, and the software that has been developed to 
analyze the data. As this software was developed, general 
purpose subroutines were also developed. Subroutines to 
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750114 THRUST (UNFIITERED 



FIGURE 3-1. THRUST TIME HISTORY 







750114 LPOP DS PR S+ 15O.01H0 
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FREQUENCY (HZ) 



manipulate and retrieve data from the data base management 
system are also described. 


3.1 Data Entry 

Initially data was provided in the form of PSD and time 
history plots. A digitizer pad was used to convert the 
plots to digitized data. Computer programs (DIGTIME and 
DIGPSD) were written to interface with the digitizer pad and 
to enter the data into the data base management system. 
Both programs function interactively with the digitizer pad 
and its operator. A complete functional description of 
these programs is provided in the following section. 

Over 300 time histories and 1200 PSDs have been digitized 
and entered into the data base management system. These 
plots were initially digitized by starting at the point 
where the trace intersected the yraxis then moving along the 
trace in straight-line segments, digitizing points between 
each segment. This resulted in data being taken at random 
time intervals for time histories and random frequency 
intervals for PSDs. The random time intervals are not 
considered to be a problem. For the PSDs, however, analysis 
efforts would be facilitated if the data was stored at 
uniform frequency intervals. The decision was made to 
convert this data to 400-line spectra to match the data 
reduction format in use at MSFC. Several methods were 
considered for making this conversion. The simplest was to 
interpolate between digitized points to obtain the necessary 
400 points at specific frequency intervals. As shown in 
Figure 3-4, this method resulted in inaccuracies in the 
magnitudes of the peaks. Several modifications of this 
method were considered, but the method that proved to be the 
most accurate utilized the following steps: 

1. Eliminate duplicated digitized points. Points 
that are within one-fourth of a bandwidth and within 10 
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percent in magnitude of each other are considered 
duplicates. 

2. For each of the 400 bandwidths in which there is 
one and only one digitized point, assign that point to 
that bandwidth. 

3. For each of the 400 bandwidths in which there are 
two or more digitized points, assign the lowest 
frequency point to the next lower bandwidth if it is 
unoccupied. Assign the second from the lowest 
frequency point within the bandwidth to that bandwidth, 
and if there is a third point, assign it to the next 
higher bandwidth if it is unoccupied. If the next 
lower bandwidth is occupied, then assign the lowest 
frequency point in the bandwidth to that bandwidth, and 
assign the highest frequency point in that bandwidth to 
the next higher bandwidth if it is unoccupied. 
Digitized points that cannot be assigned to a bandwidth 
by this method are not used. 

4. Use interpolation to obtain points for bandwidths 
that remain unoccupied. 

This conversion method was implemented in the subroutine 
CV400 . Figure 3-5 shows a comparison of digitized PSD data 
to the converted 400-line spectra. 

More recently, a method of entering data into the data base 
via magnetic tape has been developed. Time history and PSD 
data is transferred from a Norland analyzer to an HP 1000 
computer system at MSFC. The data is then put on magnetic 
tape, and the tape is read by the computer system at Wyle. 
Software was developed to read these tapes and enter the 
data directly into the data base management system. This 
method of data transfer has proven to be more accurate and 
significantly more efficient than the manual digitizing 
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• Digitized Points 
- 400 Points Interpolated 



3 - 


FIGURE 3"4. COMPARISON OF DIGITIZED POINTS TO 400-LINE 
SPECTRA OBTAINED BY INTERPOLATION 



method used initially. In this mode, spectrum interpolation 
is obviously not required. 
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Digitized Points 
400 Points Interpolated 



-5. COMPARISON OF DIGITIZED POINTS TO A00- 
SPEGTRA OBTAINED USING THE CVAOO SUBROUTINE 







3 . 2 Data Filing System 

The data base management system is designed to accommodate two types of data; time 
histories and PSD or other frequency domain data. Each time history file consists of 
two real numbers in each record. The first value in each record is a time and the 
second value is the magnitude of the quantity at that time. The data points are in 
chronological order. The file may consist of any number of data points (records), but 
most of the analysis program presently has an upper limit of 500 points. 

Rapid access to time history data is gained through a file-naming structure that is 
programmed into the data base management system. Users of the data base 
management system and its associated subroutines need not understand the file naming 
structure; however, to illustrate the flexibility of the filing system and to assist those 
attempting to expand or modify the system, the naming structure is outlined here; 

File Name; TXYYYYYY.ZZZ 
where 

T - All time history files start with the letter T. PSD files 

never start with the letter T. 

X - May be any letter or number and designates the type of 

filter used in reducing the time history data. A description 
of the filter corresponding to each letter/number is con- 
tained in a directory file named DFILT. 

YYYYYY - The 6-digit test number. 

ZZZ - A 3-digit channel number. Channel descriptors correspond- 

ing to each channel number is contained in a directory file 
named DCHN. 
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Each PSD file consists of a single real number in each record. The first N records in 
the file contain the magnitude at N evenly spaced frequencies. The last 10 records 
contain information as follows: 

Record Contents 

N+l N - the number of magnitude points. 

N+2 Composite on plots.* 

N+3 Composite computed from digitized data. 

N+4 Composite computed from data after being converted to N discrete 

points.* 

N+5 Number of digitized points that could not be used in the conversion to 

N discrete points.* 

N+6 Time slice. 

N+7 Maximum frequency. 

N+8 Spare. 

N+9 Spare. 

N+10 Spare. 

The number of discrete points, N, may be any value, but most of the analysis programs 
currently have an upper limit of 512. 


Access to PSD data files is gained through a file-naming structure similar to that used 
by the time histories. PSD files are named as follows: 

Filename: XXYYYYYY.ZZZ 
where 

XX - is a single letter followed by a single number. The letters 

A-S and U-Z and numbers 0-9 are used. The letter T is 
reserved for time histories. This letter combination desig- 
nates the maximum frequency and time slice of the PSD 
contained in that file. The time slice and maximum 
frequency corresponding to each XX is found in a directory 
file named DPSD. 

YYYYYY - The 6-digit test number. 

ZZZ - A 3-digit channel number. Channel descriptors correspond- 

ing to each channel number is contained in a directory file 
named DCHN. 


* Applies only to data obtained by digitizing plots. 
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3.3 



A total of four directory files are used by the data base management system. Three of 
these files — DFILT, DCHN, and DPSD — have already been mentioned in conjunction 
with the file-naming structure. A fourth file, named TDIR, contains the engine number, 
HPOP number, HPFP number, and the test date associated with each test number. 
Since the use of these directory files is programmed into the data base management 
system, the user need not be concerned with their use. To develop new programs, 
however, it may be necessary to know the format of the files. The format of each 
directory file is: 

TDIR - The file TDIR contains one record for each test number. Each record 
has the following format: 

Column Contents 
1-6 Test number 

8-16 Test date (YR:MO:DY) 

18-26 Engine number 

28-36 HPFP number 

38-46 HPOP number 

DCHN - The file DCHN contains a record for each channel descriptor. The 
position of the channel descriptor (i.e., record number) is the channel 
numbers. NOTE: Record must not be deleted from or inserted into 
this file. Each record has the following format: 

Column Contents 

1-40 Channel descriptor 

DFILT - The file DFILT contains a record for each filter description. Each 
record has the following format: 

Column Contents 

1 Unique letter or number for filter description 

3-15 Filter description 

DPSD - The file PSD contains one record for each PSD in the data base. The 
format of each record is as follows: 

Column Contents 

1-6 Test number 

8-10 Slice code (letter number combination used in naming PSD 
files) 

12-15 Channel number (must include high order zeros) 

16-29 Slice time 

31-44 Maximum frequency 
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3.4 S*jh£outlnea 


Subroutines that have been developed for use with the SSME data base are stored in 
FORTRAN code in a file named MSFCSUBS.FTN. A file named MSFCSUBS.OBS 
contains the assembly code for these subroutines, which are ready for use when 
establishing a program. 

To use these subroutines, simply call them in the program as required. The subroutine 
source code need not be included in the program. The following is a description of the 
subroutines included in MSFCSUBS.FTN. 

Subroutine SCHN 

Purpose: Search the channel descriptor file DCHN for a given channel descrip- 

tor and return with the corresponding channel number. 

Subroutines and functions used: None. 

Logical units used: DCHN is opened to LU 14 in SRO mode and closed before 
return. Error .messages are returned as LU 11, which must be assigned before the 
call. 

Usage: Call SCHN (ICH, CH, ISTAT) 

ICH - ASCn representation of channel descriptor in 10A4 format. ICH 
must be dimensioned in the main program to be 10 words long. 

CH - ASCII representation of the channel number in A3 format. 

ISTAT - Status on return: 0 = normal return; -1 = unable to find given 
channel descriptor — the message "CHANNEL DESCRIPTOR NOT 
FOUND" outputted to LU 11. 

Errors: 

Open Error - An error in opening DCHN results in the message "OPEN 
ERROR XXX ON DCHN" being outputted to LU 11 where XXX is the open 
error code. 
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Code Meaning 

0 No error; requested function complete. 

1 Illegal function. 

2 LU error; illegal LU. 

3 Volume error; no such volume. 

4 Name error; no such name on given volume. 

5 Size error; no room on disc for data or index block. 

6 Protect error; mismatch on protection keys. 

7 Privilege error; requested privilege may not be granted. 

8 Buffer error; no room for FCB or buffer. 

9 Assignment error; LU already assigned or off-line. 

10 File descriptor syntax error. 

12 Attempt to assign a trap generating device. 

128-255 I/O error. 


An open error results in the program being stopped. 

Close Error: An error in closing DCHN results in the message "CLOSE 

ERROR XXX ON DCHN” being outputted to LU 11 where XXX is the close 
error code. 

Code Meaning 

0 No error; requested function complete. 

2 LU error. 

9 Assignment error; LU not assigned. 

128-255 I/O error. 

A close error results in the program being stopped. 

Read Error: A read error results in the message "READ ERROR ON SCHN” 
being outputted to LU 11 and the program is stopped. 


Subroutine SPSD 

Purpose: Search the PSD directory file DPSD of a given test number, 

channel number, slice time, and maximum frequency, and return 
with the slice code. 

Subroutines and functions used: None. 

Logical units used: The file DPSD is opened to LU 13 in SRO mode and closed 
before return. Error messages as outputted to LU 11, which must be assigned 
before the call. 

Usage: Call SPSD (TEST, CH, SLICE, FREQM, ISLICE, ISTAT). 


TEST 


- ASCII representation of the 6-digit test number in A4.A2 
format. TEST must be dimensioned at least two words long. 

CH - ASCD representation of the channel number in A3 format. 

SLICE - Slice time (real variable). 

FREQM - Maximum frequency (real variable). 

ISLICE - ASCII representation of the slice code on return in A2 
format. 

ISTAT - Status on return: 0 = normal return; -1 = entry in DPSD 

matching request was not found. 

Errors: 

Open Error - An error in opening DPSD results in the message "OPEN 
ERROR XXX ON DPSD" being outputted to LU 11 where XXX is the open 
error code. 

Code Meaning 

0 No error; requested function complete. 

1 Illegal function. 

2 LU error; illegal LU. 

3 Volume error; no such volume. 

4 Name error; no such name on given volume. 

5 Size error; no room on disc for data or index block. 

6 Protect error; mismatch on protection keys. 

7 Privilege error; requested privilege may not be granted. 

8 Buffer error; no room for FCB or buffer. 

9 Assignment error; LU already assigned or off-line. 

10 File descriptor syntax error. 

12 Attempt to assign a trap generating device. 

128-255 I/O error. 

An open error results in the program being stopped. 

Close Error: An error in closing DPSD results in the message "CLOSE 

ERROR XXX ON DPSD" being outputted to LU 11 where XXX is the close 
error code. 

Code Meaning 

0 No error; requested function complete. 

2 LU error. 

9 Assignment error; LU not assigned. 

128-255 I/O error. 
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A close error results in the program being stopped. 

Read Error: A read error results in the message "READ ERROR ON SPSD" 
being outputted to LU 11 and the program is stopped. 

Subroutine DSORT 

Purpose: Sort vectors X and Y in descending order in X. 

Subroutines and Function Used: None. 

Logical Units Used: None. 

Usage: CALL DSORT (X, Y, LEN) 

X - REAL *4 vector. 

Y - REAL *4 vector. 

LEN - INTM length of vectors X and Y. 

Errors: None. 

Subroutine PEAK 

Purpose: Find the peaks in a set of data and return with the X and Y values of 

each peak. 

Subroutines Used: None. 

Logical Units Used: None. 

Usage: Call PEAK (X, Y, LEN, XP, YP, LENP). 

X - REALM X values in descending order. 

Y - REALM Y values. 

LEN - INTM length of X and Y. 

XP - REAL *4 X values of the peaks on return. XP must be 
dimensioned in calling program at least LEN/2 words long. 

YP - REALM Y values of the peaks on return. YP must be 
dimensioned in the calling program at least LEN/2 words 
long. 

LENP - INTM number of peaks found. 
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Errors: None. 


Subroutine STDIR 

Purpose: Search the test number directory TDIR for a given combination of test 

number, test date, engine number, HPFP numbr, and HPOP number. 
Any parameter(s) may be a "DON'T CARE" by calling the subroutine 
with blanks in that field. On return, all fields will be filled if search is 
successful. The subroutine may be called repeatedly to find all entries 
in TDIR that satisfy the request. 

Subroutines and Functions Called: None. 

Logical Units Used: The file TDIR is opened to LU 13 the first time the 

subroutine is called and is closed when a call results in an unsuccessful search. 
Error messages are outputted on LU 11, which must be assigned before the call. 

Usage: Call STDIR (TEST, DATE, ENG, HPFP, HPOP, ISTAT). 

TEST - ASCH representation of the 6-digit test number in A4.A2 
format. 

DATE - ASCII representation of the test date (MO:DY:YR). 

ENG - ASCII representation of the engine number in 2A4.A2 

format. 

HPFP - ASCH representation of the HPFP number in 2A4.A2 

format. 

HPOP - ASCH representation of the HPOP number in 2A4.A2 

format. 

ISTAT - Call status: 

On first call - ISTAT not equal 0. 

Call to continue search - ISTAT equal 0. 

On successful return - ISTAT equal 1. 

On unsuccessful return - ISTAT equal 0. 

On file read error - ISTAT equal 1. 
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Errors: 


Open Error - An error in opening TDIR results in the message "OPEN 
ERROR NO XXX" being outputted to LU 11 where XXX is the open error 
code. 

Code Meaning 

0 
1 
2 

3 

4 

5 

6 

7 

8 

9 

10 
12 

128-255 

An open error results in the program being stopped. 

Subroutine RDTH 

Purpose: Read the time history file specified by test number, filter code, and 

channel number, and return with the contents. 

Subroutines and Functions Used: None. 

Logical Units Used: The specified time history file is opened to LU 12 and closed 
before return error messages are outputted to LU 11, which must be assigned 
before the call. 

Usage: Call RDTH (TEST, CH, IFLT, TIME, MAG, LEN, 1ST AT). 

TEST - ASCII representation of the test number in A4.A2 format. 

CH - ASCH representation of the channel number in A3 format. 

IFLT - ASCH representation of the filter code in A1 format. 

TIME - REAL *4 time values on return. 

MAG - REAL *4 magnitude values on return 

LEN - Dimension of vectors TIME and MAG on CALL number of 
data points read on return. 


No error; requested function complete. 

Illegal function. 

LU error; illegal LU. 

Volume error; no such volume. 

Name error; no such name on given volume. 

Size error; no room on disc for data or index block. 
Protect error; mismatch on protection keys. 

Privilege error; requested privilege may not be granted. 
Buffer error; no room for FCB or buffer. 

Assignment error; LU already assigned or off-line. 

File descriptor syntax error. 

Attempt to assign a trap generating device. 

I/O error. 
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1ST AT - 0 = Normal return 

-1 = Number of data points greater the LEN at CALL 

1 = Illegal function 

2 = LU error; illegal LU 

3 = Volume error; no such volume 

4 = Name error; no such mane on given volume 

5 = Size error; no room on disc for data or index block 

6 = Protect error; mismatch on protection keys 

7 = Privilege error; requested privilege may not be granted 

8 = Buffer error; no room for FCB or buffer 

9 = Assignment error; LU already assigned or off-line 

11 = File descriptor syntax error 

12 = Attempt to assign a trap generating device 
128-255 = I/O error 

Errors: 

Read Error - A read error results in the message "Read Error on File XXX" 
being outputted to LU 11 where XXX is the file name being read. 

Subroutine ERROR 

Purpose: Output an error message on LV11 in the form: 

Error message ISTAT = 10 

Subroutines and Functions Used: None. 

Logical Units Used: The message is outputted to LU 11, which must be assigned 
before the call. 

Usage: Call ERROR (MESG, ISTAT). 

MESG - May be a 1X5 array with ASCII representation of the 
message or the message may be placed in the call CALL 
ERROR ("MESSAGE GOES HERE: ISTAT) the message 
must be 20 characters long including blanks ISTAT - INT*4 
variable or constant to be outputted. Example: ISTAT = 
10; IF (ISTAT. NE.0) CALL ERROR ("THIS IS A SAMPLE ' 
ISTAT) will result in the following message being outputted 
to LU 11: THIS IS A SAMPLE ISTAT = 10. 
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Errors: None. 


Subroutine PROMPT 

Purpose: Output a message to LU11. 

Usage: Call PROMPT(TEXT) 

TEST - Message is in the form "EXAMPLE MESSAGE!" the mes- 
sage may be any length but must end with a !. Example: 
Call PROMPT ("EXAMPLE MESSAGE!") result in the fol- 
lowing being outputted to LU 11: EXAMPLE MESSAGE!. 

Errors: None. 

Subroutine INTERP 

Purpose: Interpolate between data points. 

Usage CALL INTERP (X, Y, XX YY, SLOPE, LEN). 

X - REAL *4 array containing x data points. 

Y - REAL *4 array containing y data points. 

XX - REAL *4 scalar the x value for which a y value is to be 

computed. 

YY - REAL *4 scalar y value found to correspond to given x 
value. 

SLOPE - REAL *4 scalar slope at XX. 

LEN - The length of the vectors X and y. 

Errors: None. 

Subroutine RDPSD 

Purpose: Read the contents of a PSD data file and put frequency values in X 

and magnitudes in Y. 

Subroutines and Function Called: None. 
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Logical Units Used: The PSD data file is opened to LU 12 and closed before the 
return usage: Call RDPSD (TEST, CH, ISLICE, X, Y, COMP, LEN, ISTAT) 

TEST - Test number in A4-A2 format. 

CH - Channel number in A3 format. 

ISLICE - Slice code in A2 format. 

X - REAL *4 any containing frequency value on return must be 

dimensioned at least LEN words long in calling program. 

Y - REAL *4 array containing magnitude values on return must 

be dimensioned at least LEN+10 words long in calling 
program. On return: 

Y(LEN+1) = LEN 

Y(LEN+2) = composite on plots 

Y(LEN+3) = composite computed from digitized data 

Y(LEN+4) = composite computed from discrete spectra 

data 

Y(LEN+5) = number digitized points not assigned during 
conversion to discrete spectra 
Y(LEN+6) = time slice 
Y(LEN+8)-Y(LEN+10) = spare 

LEN - Dimension of X at call and number of data points read on 
return. 

ISTAT - Call status on return 

-625 to -755 I/O error at close 

-509 Assignment error; LU not assigned at close 

-502 LU error at close 

-3 LEN check failed 

-2 Data exceed dimension of X 

-1 Read error 

0 Normal return 

1 Illegal function 

2 LU error; illegal LU 

3 Volume error; no such volume 
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4 Name error; no such name on given volume 

5 Size error; no room on disc for data or index block 

6 Protect error; mismatch on protection keys 

7 Privilege error; requested privilege may not be granted 

8 Buffer error; no room for FCB or buffer 

9 Assignment error; LU already assigned or off-line 

11 File descriptor syntax error 

12 Attempt to assign a trap generating device 
128-255 I/O error 


Errors: None. 


Subroutine SPSD1 

Purpose: Search the PSD directory file for a given test number and channel 

number, and return with the time slice maximum frequency and slice 
code. The subroutine can be called repeatedly to find all entries that 
satisfy the request. 

Subroutines and Function Called: None. 

Logical Unit Used: The DPSD file is opened to LU 13 and remains open until a 
search is unsuccessful. Error messages are outputted on LU 11, which must be 
assigned before the call. 


Subroutine CV400 

Purpose: Converts digitized PSD data taken at nonuniform frequency intervals 

to 400 data points at constant frequency intervals. 
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3.5 Programs 


Both data base management and data analysis programs have been written and used 
with the SSME data. The data base management program includes programs for 
entering new data, for summarizing the contents of the data base, and for removing 
erroneous or incorrectly identified data. Data analysis programs include programs for 
plotting data and computing statistcal data summaries. The following is a brief 
summary of the programs that have been completed to date. 

Program name: DIGTIME 

Purpose: Used in digitizing time history plots. 

Method: The program asks the operator to enter the filter description, test 

number, and channel descriptor for the plot to be digitized. It then 
checks the TDIR directory file to see if the test number has already 
been entered. If the entry is found, the test date, engine number, 
HPOP number, and HPFP number is displayed for verification. If no 
entry is found, the operator is asked to enter the test date, engine 
number, HPOP number, and HPFP number, and this information is 
entered into TDIR. 

The DCHN directory file is then searched for the channel descriptor 
and to obtain the corresponding channel number. If the channel 
descriptor is not found, then the program asks if it is a new channel 
descriptor. If the operator responds with "yes", then an entry is made 
in DCHN assigning to that channel descriptor a channel number. If the 
response is "no" then the operator is given the opportunity to reenter 
the channel descriptor. 

The program then asks the operator to enter the X and Y values at the 
origin of the plot to be digitized and the values of one other point on 
the x-axis and one on the y-axis. It then instructs the operator to 
digitize these points using the digitizer pad. From this information 
the program determines the orientation of the plot on the digitizer pad 
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and the scaling parameter needed to convert the output of the 
digitizer pad into engineering units. The operator digitizes the plot in 
as many straight-line segments as is necessary to produce an accurate 
digital image. When the digitizing is complete, the program creates a 
file, using the naming structure described previously, and enters the 
data in that file, completing the process of entering a time history 
into the data base. 

Program name: DIGPSD 

Purpose: Used in digitizing PSD plots. 

Method: The DIGPSD program starts by asking the operator for. the test number 

and channel descriptor — much like DIG TIME — and, if necessary, makes 
the same search of entries in the TDIR and DCHN directory files. In 
addition, it asks for the time slice and maximum frequency of the PSD 
to be digitized. With this information, it searches the DPSD directory 
file to find a unique slice code for the plot to be digitized. The slice 
code, test number, and channel number are used to name the data file 
to be created. The program asks the operator to digitize three points 
on the axes of the plot from which it determines the plot orientation 
and scaling factors. The program assumes that the magnitude axis is 
1'Qgorlthm and the frequency axis is linear. After the plot is digitized, 
the program converts the data to 400 magnitude points at uniform 
frequency intervals. These points are stored in a data file named using 
the file naming structure described previously. The program asks the 
operator to enter the composite value. It also computes the composite 
from the digitized data and the data after being converted to 400 
points. All three composite values are entered into the data file, and 
the computed values are displayed so the operator can check the 
digitization. An entry is then made in the DPSD directory file, 
completing the data entry process. 

Program name: REMTH 

Purpose: Remove time history data from the data base. 
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.Method: 


Program 

Purpose: 

Method: 


Program 

Purpose: 

Method: 


Program 

Purpose: 


The program asks for the test number, filter description, and channel 
descriptor of the time history to be removed. It then reads the 
corresponding time history file and displays the minimum and maxi- 
mum values of magnitude and time as well as the number of data 
points in the file. The operator is then asked, by the program, if this 
is the file to be removed. If the operator's response is yes, then the 
file is removed; otherwise it remains unchanged. 

name: REMPSD 

Remove PSD data from the data base. 

The program asks the operator to enter the test number, channel 
descriptor, time slice, and maximum frequency of the PSD data to be 
removed. It then searches the DPSD directory file for the given 
combination of test number, channel descriptor, time slice, and 
maximum frequency to obtain the slice code. The data file is then 
read. The minimum and maximum magnitudes and composite value is 
displayed. The operator is then asked if this is the PSD data to be 
removed. A response of yes by the operator results in the data file 
being deleted and the corresponding entry in the DPSD directory file 
being removed. Any other response results in no action being taken. 

name: LSTTH 

List the test number and channel descriptor of all time history data in 
the data base. 

The program uses a system command to obtain a list of all the time 
history file names from which the test numbers and channel numbers 
are obtained. The DCHN directory file is read to obtain the channel 
descriptor corresponding to each channel number. The test number 
and channel descriptor for each data file is then printed. 

o 

name: LSTPSD 

List the test number, time slice, maximum frequency, and channel 
descriptor of each PSD in the data base. 
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Method: 


The DPSD directory file is read to obtain the test number, time slice, 
maximum frequency, and channel number of each PSD in the data 
base. The DCHN directory file is used to find the channel descriptor 
that corresponds to each channel number. The test number, time 
slice, maximum frequency, and channel descriptor is then printed. 

Program name: PLOTTH 

Purpose: Plot time history data. 

Method: The program asks for the test number filter description and channel 

descriptor. The channel number corresponding to the given channel 
descriptor is found by searching the DCHN directory file. The data 
file is read and the data is plotted on a Tektronix graphics terminal 
from which a hard copy may be produced. The plot includes scaled and 
labeled axes and a title. The title consists of the test number, channel 
descriptor, and filter description. 

Program name: PLOTPSD 

Purpose: Plot PSD data. 

Method: The program asks the operator to enter the test number, channel 

descriptor, time slice, and maximum frequency of the PSD data to be 
plotted. The program searches the DCHN directory file to find the 
channel number and the DPSD directory file for the slice code. The 
corresponding data file is read and the data is plotted on a Tektronix 
graphics terminal. The axes are scaled and labeled, and a title 
consisting of the test number, channel descriptor, and slice time is 
included on the plot. 

Program name: STATPSD 

Purpose: For a set of PSDs specified by engine number, HPFP number, HPOP 

number (any or all may be "don't care"), maximum frequency, channel 
descriptor, and thrust range, calculates the following as a function of 
frequency. 

- Mean 
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- Standard deviation 

- rms 

- Third central moment 

- Fourth central moment 

- Minimum value 

- Maximum value 

- Range 

Method: The program asks the operator to enter the engine number, HPOP 

number, and HPFP number of the data to be considered. Entries that 
are left blank are assumed to be "don't care." The program searches 
the TDIR directory file to find the test numbers that correspond to the 
given combination of engine number, HPOP number, and HPFP num- 
ber. The DCHN directory file is searched to find the channel number 
corresponding to the channel descriptor. The DPSD file is searched to 
find slice codes and time slices that correspond to the given channel 
number, maximum frequency, and each test number. For each entry 
found, the corresponding thrust time history is read to find the thrust 
at that time slice. If the thrust is within the specified range, then the 
PSD data file is read and included in the statistical calculations. Once 
all the PSDs that correspond to the given set of requirements have 
been found, the data is outputted to a file in a format to be used by 
other programs. This output file also includes information needed to 
title and label the axis of plots. 

Program name: EXPSD 

Purpose: For each PSD as specified by engine number, HPFP number, HPOP 

number (any or all may be "don't care"), maximum frequency, and 
channel descriptor, finds the magnitude and frequency of the N largest 
peaks, where N is specified by the operator. Also, finds the thrust 
that corresponds with each PSD considered. The output of this 
program can be used to plot the magnitude or frequency of a given 
peak as a function of thrust. 
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Method: PSD files that satisfy the specification are found by using the method 

described for STATPSD. The peaks in each PSD are found by 

examining the slope of the PSD as a function of frequency. When the 
slope changes from positive to negative, that point is stored as a peak. 
When all the peaks in a PSD are found, they are sorted by magnitude, 
and the N largest peaks are reported to an output file together with 
the thrust corresponding to that PSD. This process is repeated for 
each PSD that satifies the given specification. 

Program name: PLOT 

Purpose: Plot the contents of a file generated by analysis programs such as 

STATPSD and EXPSD. 

Method: The data file to be plotted is specified by name when the program is 

started. From this file, the program reads a description of each 

channel of data that the file contains. These descriptions are 
displayed and the operator is asked to enter a number corresponding to 
the channel to be plotted on the x-axis and for the channel to be 
plotted on the y-axis. The operator is asked to specify linear or 
logarithmic scaling for each axis. The choice of line plot or dots is 
given. If line plot is specified, the option of sorting the data by x- 
values is given. The plots are titled and axes are labeled using the 
information provided by the data file. 

Program name: LOADTHST 

Purpose: Enter thrust profile or rms time history data into the data base via 9- 

track magnetic tape. 

Method: The program reads the test number, engine number, HPOP number, 

HPFP number, and test data from the tape. The TDIR file is checked 
for the test number. If the test number is not found, the test number, 
HPOP number and HPFP number are added to the file. The program 
then reads the time history from the tape and looks to see whether the 
information is already in the data base. If it is not in the data base, 
then the history read from the tape is added to the data base. 
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Program name: LOADPSD 

Purpose: Enter PSD data into the data base. 

Method: The program reads the test number channel descriptor and time slice 

from the tape. It then checks the TDIR file to see if that test 
number is present. If not found, the program prints an error 
message and proceeds to the next data set on the tape. The 
program then checks the DPSD file for a data file with the same 
test number, channel descriptor, and time slice. If one is 
found, the program proceeds to the next data set on the tape; 
if no match is found, a slice code is obtained and the PSD data 
is transferred from tape to the data base management system. 

A flow diagram of the data base management and analysis software is 
illustrated in Figure 3-6. 

A computer listing of the programs and subroutines discussed in this 
section is provided on the following pages. This code is operational 
on PC/compatible systems. The routines are also available on disc for 
efficient transfer to MSFC computers. 
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09:32:46 

1 7 Microsoft F0RTRAN77 V3. 20 02/84 

*DEBUG 

C MAIN PROGRAM TO CALCULATE STATISTICS ON PSD’ S 

C THE PROGRAM SEARCHES THE DATA BASE FDR A GIVEN COMBINATION 

C OF ENGINE NO., HPFP NO., HPOP NO. (ANY OR ALL MAY BE "DON’T CARE") 

C CHANNEL DESCRIPTOR AND THRUST RANGE. 

C 

C THE MEAN, STANDARD DEV I AT ION, TH I RD CENTRAL MOMENT, FOURTH CENTRAL 

C MOMENT, RMS, MINIMUM VALUE, MAX I MUN VALUE AND RANGE ARE CALCULATED 

C 

C OUTPUT DATA IS WRITTEN ON LU 7 (120 CHARACTORS) 

C THE GIVEN ENGINE NUMBER, HPFP NUMBER AND HPOP NUMBER, TEST NUMBERS 

C PSD’S AND QUANTITY OF PSD FROM EACH TEST USED ARE OUT PUT ON LU 5 

C 
C 
C 

CHARACTER TEST*6, DATE*8, ENG*10, HPFP*10, HP0P*10 
CHARACTER TESTX*6, DATEX*8, ENGX*10, HPFPX*10, HPOPX* 10 
CHARACTER ICH*£0, CH*3 , 

REAL TIME (500) , MAG (500) , X ( 500) , Y (500) 

REAL XXFOUR (200) 

REAL XXTOT (200) , XXMIN(200) , XXMAX (200) , XXSQR(EOO) , XXCUB(200) 

OPEN (7, F I LE=’ DATA. PRT’ ) 

OPEN (5, FILE=’ STATS. DAT’ ) 

C 

C INITIALIZE VARIABLES 

C 


27 

1 28 

1 29 

1 30 

1 31 

1 32 

1 33 5 

34 

35 

36 20 

37 22 

38 30 

39 

40 1 0 

41 C 

42 C 

43 C 

44 1 

45 40 

46 

47 


DO 5 J=l, 200 
XXTOT (J) =0.0 
X XM I N ( J ) = 1 . 0E32 
XXMAX (J)=-l. 0E32 
XXCUB ( J) =0. 0 
XXFOUR (J) =0. O 
X XSQR ( J) =0. 0 
I ST AT=0 
I ST ATT =0 
FORMAT (A10) 
FORMAT (’ ’A10) 

FORMAT ( A20 ) 
TEST=’ ’ 

DATE = ’ » 


ASK FOR ENGINE NO. , HPFP NO. , HPOP NO. 
WRITE (*,40) 

FORMAT (’ ENTER ENGINE NO.’) 

READ (*,20) ENG 
IF (ENG. EQ. ’ STOP’ ) STOP 


48 WRITE (*,50) 

49 50 FORMAT (’ ENTER HPFP NO. ’) 

50 READ <*, 20) HPFP 

51 WRITE (*,60) 

52 60 FORMAT (’ ENTER HPOP NO. ’ ) 

53 READ (*,20) HPOP 

54 C 

55 C ASK FOR CHANNEL DESCRIPTOR AND FIND CHANNEL NO. 

56 C 3-32 



Zj / 


WKi 1 C. <■*, ' fclN 1 tr< UI-IHININtli- ' 




5 a 


READ <*, 30) ICH 




59 


CALL SCHIMC ICH, CH, ISTAT) 
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60 


IF C ISTAT. EQ. -1) GO TO 1 




61 


IF ( ISTAT. NE. 0) CALL ERROR ( ’ IN ICHN ISTAT) 




62 


CALL DPSDTEMP(TESTX, CH, ISTAT) 




63 

C 





64 

c 

ASK FOR MAXIMUM FREQUENCY 




65 

c 





66 


WRITE <*,*)’ ENTER MAX. FREQUENCY!’ 




67 


READ (*, 21 ) FREQM 




66 


AFREQS-FREQM/400. 




69 

ro 

FORMAT (G13. 6) 




70 


FREQL=FREQM 




71 

c 





72 

c 

ASK FOR THRUST RANGE 




73 

c 





74 


WRITE (*,*)’ ENTER LOWER THRUST LIMIT!’ 




75 


READ <*,£!) THRL 




76 


WRITE**, *)’ ENTER UPPER THRUST LIMIT!’ 




77 


READ (*, 21 ) THRH 




78 

c 





79 

c 

WRITE NUMBER OF OUT PUT COLUMNS TO LU 7 (THIS CASE 9) 




60 

c 





61 


WRITE (7, 64) 




62 

64 

FORMAT (’ 0009’) 




83 

C 





84 

C 

WRITE HEADER 




85 

C 





86 


WRITE (7, 65) ICH, THRL, THRH 




87 

65 

FORMAT (’ ’ A10, 5X, ’ THRUST RANGE = ’ F5. 1,’ - ’ F5. 1) 




88 

c 





89 

c 

WRITE OUTPUT COLUMN ID’S 




90 

c 





91 


WRITE (7, 66) 




92 

66 

FORMAT ( ’ FREQUENCY HZ’/’ MEAN P**2/HZ’ 




93 


1/’ STANDARD DIVIATION P**2/HZ’ /’ THIRD CENTRAL MOMENT 

P**6/HZ’ 



94 


4/’ FOURTH CENTRAL MOMENT P**8/HZ’ 




95 


£/’ RMS P**2/HZ’/’ MINIMUM VALUE P**2/HZ’ 




96 


3/’ MAXIMUM VALUE P**2/HZ’ /’ RANGE P**2/HZ’ ) 




97 

C 





98 

c 

WRITE ENGINE NO., HPFP NO., HPOP NO. COMBO TO LU 5 




99 

c 





100 


WRITE (5, 67) ENG 




101 


WRITE (5, 68) HPFP 




102 


WRITE (5, 69) HPOP 




103 

67 

FORMAT (’ ENGINE NO. = ’ A10) 




104 

68 

FORMAT (’ HPFP NO. = ’ A10) 




105 

69 

FORMAT <’ HPOP NO. = ’A10 //’ TEST NO. QTY. ’ ) 




106 


JK=0 




107 

C 





o 

CD 

C 

SEARCH FOR TEST NO. FOR GIVEN COMBO 




109 

C 





1 10 

70 

CONTINUE 




111 


TES7X=TES7 




i 12 

71 

DATEX=DATE 




113 


ENGX=ENG 




114 


HPOPX=HPOP 




115 

76 

WPFPX=HPFP 
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IF (ISTATT. LT. 0) GO TO 930 
IF ( ISTATT. EQ. 0) GO TO 300 


n Hr PA, DHUKA, 
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Line# 
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1 19 



120 



121 

C 


122 

C 


123 
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157 



158 

c 


159 

c 


160 

c 


161 
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7 Microsoft F0RTRAN77 03.20 02/84 

NC0UT=0 

LENT — 1 00 

READ THRUST DATA 

CALL RDTH (TESTX, ’ 001 ’,’1 ’ , T IME, MAG, LENT, I STAT) 

IF ( I ST AT. EQ. 4) GO TO 1100 

IF ( ISTAT. NE. 0) CALL ERROR (’IN RDTH ’ , ISTAT) 

I STAT 1=1 

SEARCH FOR PSD DATA 

CALL SPSD1 (TESTX, CH, SLICE, FREQ, ISLICE, ISTAT1, THRUST) 

IF FILE FOUND GO TO 110 
IF ( ISTAT 1 . NE. —1 ) GO TO 110 

IF NO MORE FOUND WRITE TEST NO. AND QUANTITY OF PSD’ S FOUND 
FOR THAT TEST NO. . . GET NEW TEST NO. 

IF (NCOUT. EQ. 0) GO TO 70 
WRITE ( *, 105) TESTX, NCOUT 
WRI TE < 5, 105) TESTX, NCOUT 
FORMAT < ’ ’ A6, 2X, 13) 

GO TO 70 

CONTINUE 

IF OUT SIDE THRUST RANGE GET NEW PSD 

IF ( (THRUST. LT. THRL) . OR. (THRUST. GT. THRH) ) GO TO 100 

IF PSD DOES NOT COVER FREQUENCY RANGE GET NEW PSD 

IF (FREQ. NE. FREQL) GO TO 100 

NCOUT =NCQUT + 1 

LEN=500 


READ PSD DATA FILE 

CALL RDPSD (TESTX, CH, ISLICE, X, Y, COM, LEN, ISTAT) 

IF ( ISTAT. NE.O) CALL ERROR (’IN RDPSD ’, ISTAT) 

JK= JK+ 1 

DO 200 J= 1,400 
FREQZ= ( J-l ) *AFREQS 

ACCUMULATE DATA FOR 400 FREQUENCIES BETWEEN O AND MAX FREQUENCY 

XXTOT ( J) =XXTOT ( J) +Y < J) 

XXMIN(J) = AM INI ( XXMIN ( J) , Y ( J) ) 

XXMAX ( J) =AMAX1 (XXMAX (J) , Y(J) ) 

XXSQR ( J ) =X XSQR ( J) +Y( J) **2 
XXFOUR ( J) =XXFOUR ( J) +Y < J) **4 
XXCUB ( J) =XXCUB ( J) +Y ( J ) ***3 


i 


X i/O cuu 

UUIN 1 1 »NUC 

176 

GO TO 100 

177 300 

CONTINUE 

D Line# 1 
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178 

IF < JK. GT. 0) GO TO 301 

179 

WRITE <*,*)’ NONE FOUND!’ 

ISO 

STOP 

181 C 

188 C 

CALCULATE AND OUTPUT STATISTICS j 

1 183 C 

184 301 

DO 400 J=l,400 

1 185 

FREQS= ( J-l ) *AFREQS 

1 186 

XMEAN=XXTOT ( J) /JK 

1 187 

XM£=X XSQR ( J) /JK 

1 188 

VAR=XM£-XMEAN**£ 

1 189 

STD=SQRT ( ABS (VAR) ) 

1 190 

RMS=SQRT (XM2) 

1 191 

RANGE=X XMAX ( J)-XXMIN(J) 

1 198 

XM3=X XCUB ( J ) /JK 

1 193 

XM4=XXF0UR < J) /JK 

1 194 

XU3— XM3— 3. *XMEAN*XM2— 2- *XMEAN**3 

1 195 

XU4=XM4— 4. *XM3*XMEAN+6. *XM2*XMEAN**2-3. *XMEAN**4 

1 196 400 

WRITE (7, 310) FREQS, XMEAN, VAR, STD, XU3, RMS, XXMIN < J) , XXMAX < J) 

1 197 1 

, RANGE 

198 310 

FORMAT (’ ’ 9G13. 6) 

199 

WRITE ( 5, 320) JK 

£00 3£0 

FORMAT (’ TOTAL ’13) 

201 

STOP 

202 990 

WRITE (*,991) 

£03 991 

FORMAT <’ FILE READ ERROR’) 

£04 

STOP 

205 1100 

WRITE ( *, 1110) TESTX 

206 1110 

FORMAT (’ THRUST DATA FOR TEST NO. ’ A4, A£ 

207 1 

’ NOT IN DATA BASE’ ) 

£08 

GO TO 70 

£09 

END j 

Name Type 

Offset P Class 

ABS 

INTRINSIC 

AFREOS REAL 

1 3000 

AMAX1 

INTRINSIC 

AMIN! 

INTRINSIC 

CH CHAR*3 

12986 

COM REAL 

13640 

DATE CHAR*S 

12859 

DATEX CHAR*8 

13556 

ENG CHAR* 10 

12891 

ENGX CHAR* 10 

13564 

FREQ REAL 

13610 

FREQL REAL 

13012 

FREQM REAL 

18996 

FREQS REAL 

13648 

FREQZ REAL 

13644 

h !D-p CHAR* 10 

12924 

HPFPX CHAR *10 

13584 

h PO ^ CHAR* 1 0 

12956 

HPOPX CHAR* 10 

13574 

ICh CHAR *20 

12966 

ISuICE INTEG£R*4 13614 




.Ltoim l railed* *+ 1 c'dcu 

ISTAT1 INTEGER*4 13602 

ISTAT7 INTEGER*4 12824 


D Line# 1 7 

J 

INTEGER*4 

JK 

INTEGER*4 

LEN 

INTEGER*4 

LENT 

INTEGER*4 

MAG 

REAL 

NCOUT 

INTEGER*4 

RANGE 

REAL 

RMS 

REAL 

SLICE 

SORT 

REAL 

STD 

REAL 

TEST 

CHAR*6 

TESTX 

CHAR*6 

THRH 

REAL 

THRL 

REAL 

THRUST 

REAL 

TIME 

REAL 

VAR 

REAL 

X 

REAL 

XM2 

REAL 

XM3 

REAL 

XM4 

REAL 

XMEAN 

REAL 

XU3 

REAL 

XU4 

REAL 

XXCUB 

REAL 

XXFOUR 

REAL 

XXMAX 

REAL 

XXMIN 

REAL 

XXSQR 

REAL 

XXTOT 

REAL 

Y 

REAL 


12Q16 
13552 
1 3636 
13598 
10616 
13594 
13672 
13668 
13606 

INTRINSIC 

13664 
12853 
12989 
13020 
13016 
13618 
16 
13660 
2016 
13656 
13676 
1 3680 
1 3652 
13684 
13688 
10016 
6016 
9216 
7616 
8416 
6816 
4016 
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210 
21 1 
212 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 


24 



226 

227 

228 

229 

230 


C 

C*PR0G SCHN 

C PURPOSE: SEARCH THE DESCRIPTOR FILE FOR A GIVEN DESCRIPTOR 

C AND RETURN THE CHANNEL NUMBER 

C 

SUBROUT I NE SCHN < I CH, CH, I STAT ) 

C 

C ICH - DISCRIPTOR (20 CHARACTERS) 

C CH CHANNEL NO- (RETURNED) 

C ISTAT - RETURN STATUS 

C O-NORMAL 

C — 1— DISCRI PTOR NOT FOUND 

CHARACTER ICHY*4 

CHARACTER FD14*4, CH*3, ICH*20, ICHX*20, DMY*1 
DATA FD14/’ DCHN 1 / 

OPEN (14, FIL£=FD 1 4, I OSTAT= I STAT ) 

IF ( ISTAT- NE. O) GO TO 2040 
I =0 

3000 1 = 1 + 1 

READ (14, 703, END=3500, ERR= 1000) ICHX 
703 FORMAT (A20) 

IF ( ICHX. NE. ICH) GO T 0 3000 


£34 


UL.U£>t ^ i *+, iuaiHi-icrtrt; 

IF ( I ERR- NE. O) GO TO £060 
1 = 1+1000 


D Line# 

1 

7 Microsoft 

£35 


WRITE (ICHY, 3030) I 

£36 

3030 

FORMAT (14) 

£37 


GO TO 3600 

£38 

3500 

WRITE (•*, 3501) 

239 

3501 

FORMAT (’ CHANNEL DESCRIPTOR NOT FOUND’) 

£40 


I STAT =— 1 

241 


CLOSE (14, IOSTAT=IERR) 

242 


IF ( IERR. NE. 0) GO TO £060 

£43 


RETURN 

244 

3600 

READ < ICHY, 701 ) DMY, CH 

£45 

701 

FORMAT (Al, A3) 

246 


WR I TE ( *, 701 ) DMY, CH 

247 


RETURN 

£48 

1000 

WRITE ( •*, 1110) 

249 

1110 

FORMAT ( ’ READ ERROR ON SCHN’ ) 

£50 


STOP 

£5 1 

2040 

WRITE <*, £050) ISTAT 

252 

£050 

FORMAT <’ OPEN ERROR ’13’ ON DCHN’ ) 

253 


STOP 

254 

£060 

WR I TE ( *, £070) IERR 

£55 

£070 

FORMAT (’ CLOSE ERROR ’13’ ON DCHN’) 

£56 


STOP 

£57 


END 
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Name 

Type 

CH 

CHAR*3 

DMY 

CHAR*1 

FD 1 4 

CHAR** 

‘I 

INTEGER** 

ICH 

CHAR*20 

I CHX 

CHAR*20 

ICHY 

CHAR** 

IERR 

INTEGER** 

ISTAT 

INTEGER** 


Offset P Class 

4 * 

13904 
1 3Q£4 
13828 
O * 

13832 

13862 

13858 

8 * 


£58 C*PROG SPSD 

253 C PURPOSE: SEARCH THE PSD DIRECTORY FOR A GIVEN TEST NO. 

260 C CHANNEL NO. , TIME SLICE, AND MAX FREQ.- — RETURN 

261 C WITH THE SLICE CODE 

262 C NOT VERIFIED ON IBM 12/85 
£63 C 


£64 

£65 

£66 

£67 

£68 

£69 

£70 6£0 
£71 

c. / £ SjO 
£73 

£74 840 

£75 

£76 


SUBROUTINE SPSD < TEST, CH, SL I CE, FREQM, ISLICE, ISTAT) 
CHARACTER TEST*6, TEST1*6 
CHARACTER FD13*4 
DATA FD1 3 /’ DPSD’ / 

OPEN (13, FI LE=FD 13, I OSTAT= I STAT ) 

IF ( ISTAT. NE. 0) GO TO 1180 

READ (13, 830, 2ND=S60, ERR=1000) TEST 1 , ISLICE, CHX, SLICE1, 
1 FREQ 

FORMAT <A6, IX, A£, IX, A3, £(1X, G13. 6) ) 

IF (TEST. NE. TEST1 ) GO TO 820 
CONTINUE 

IF (CHX. NE. CH) GO TO 8£0 
IF (SLICE. NE. SLICED GO TO 8£0 


£78 850 
£79 


11“ V f Ktlb!. iMt. ntzUVU U3U lU tfc.U 

CLOSE (13, IOSTAT=IERR) 

IF < I ERR. NE. O) GO TO £000 


Line# 

£80 

£81 

£ 8 £ 

£83 

£84 

£85 

£86 

£87 

£88 

£89 

£90 

£91 

£9£ 


100O 

1110 

1180 

1190 

£000 

£010 


7 

RETURN 
ISTAT =— 1 
GO TO 850 
WRITE <*, 1110) 

FORMAT ( 1 READ ERROR IN £ 
STOP 

WRITE <*, 1190) ISTAT 
FORMAT <’ OPEN ERROR NO. 
STOP 

WRITE (*, £010) I ERR 
FORMAT ( 1 CLOSE ERROR NO. 
STOP 
END 
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IN SPSD’ ) 


’ 13, * 


FILE DSPD’ ) 


’ I 3’ ON FILE DPSD’ ) 


Name Type 


Offset P Class 


CH 

REAL 

4 

* 

CHX 

REAL 

14006 


FD1 3 

CHAR*4 

13996 


FREQ 

REAL 

14014 


FREQM 

REAL 

1 £ 

* 

I ERR 

integer*^ 

14046 


I SLICE 

INTEGER*^ 

16 

* 

ISTAT 

INTEGERS 

£0 

* 

SLICE 

REAL 

8 

* 

SLICE1 

REAL 

14010 


TEST 

CHAR*6 

O 

* 

TEST 1 

CHARGE 

14000 




£94 

C 



£95 

C 

PURPOSE: SORT X AND Y IN 


£96 

c 

NOT VERIFIED ON IBM l£/85 


£97 

c 



£98 


SUBROUTINE DSORT ( X, Y, LEN) 


£99 


DIMENSION X(1),Y(1) 


300 


LM1 =LEN-1 


301 


DO 10 1=1, LM1 

1 

30£ 


LEND=LEN— I 

1 

303 


DO 10 J=1 , LEND 

£ 

304 


JP1=J+1 

£ 

305 


IF ( X ( J) . GT. X ( JP1 ) ) GO TO 10 

£ 

306 


XTEMP=X (JP1) 

£ 

307 


YTEMP=Y< JP1) 

c. 

308 


X ( JP1 ) =X ( J) 

£ 

309 


Y ( JP 1 ) =Y < J ) 

£ 

310 


X ( J) =XTEMP 

p 

311 


Y(J)=YTEMP 

£ 

31 £ 

10 

CONTINUE 


313 


END 


TO DECENDING ORDER 


Offset 


C 1 ass 


IN i EG£R*4 
TNTEGER*4 


14156 

14168 


JW1 J.IM I fc.lafc.H'*'* iH-i/O 

LEN INTEGER*4 8 * 

LEND INTEGERS 14164 
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LM1 INTEGER*4 14152 

X REAL 0 * 

XTEMP REAL 14180 

Y REAL 4 * 

YTEMP REAL 14184 


314 C*PR0G PEAK 

315 C PURPOSE: SEARCH PSD DATA TO FIND PEAKS 

316 C NOT VERIFIED ON IBM 12/85 

317 C 

318 SUBROUTINE PEAK ( X, Y, LEN, XP, YP, LENP) 


319 

C 



320 

C 

X 

PSD FREQ POINTS 

321 

c 

V 

PSD MAG. POINTS 

322 

c 

LEN - 

NUMBER OF DATA POINTS 

323 

c 

XP 

- FREQUENCY OF PEAKS 

324 

c 

YP - 

MAGNITUDE OF PEAKS 

325 

c 

LENP 

- NUMBER OF PEAKS 

326 

c 




327 DIMENSION X ( 1 ) , Y ( 1 ) , XP ( 1 > , YP ( 1 > 

328 LENM1=LEN— 1 

329 LENP=0 

330 DO 10 J=2, LENM1 

1 ' 331 JP1=J+1 

1 332 JM1=J— 1 

1 333 IF(. NOT. (Y( J)-Y ( JM1 ) . GT. O. O) . AND. (Y ( J)-Y( JP1 ) . LT. O. O) ) GO TO 10 

1 334 LENP=LENP+1 

1 335 XP (LENP) =X < J) 

1 336 YP (LENP) =Y ( J> 

1 337 10 CONTINUE 

338 RETURN 

339 END 


Name 

Type 

Offset 

P 

J 

INTEGER*4 

14192 


JM1 

INTEGER*4 

14204 


JP1 

INTEGERS 

14200 


LEN 

INTEGER*4 

a 

* 

LENM1 

INTEGER*4 

14188 


LENP 

INTEGERS 

20 

* 

X 

REAL 

0 

* 

XP 

REAL 

12 

* 

Y 

REAL 

4 

* 

YP 

REAL 

16 

* 


340 C 

341 C3PR0G STDIR 

342 C SUBROUTINE STDIR 

343 C 

344 C PURPOSE: SEARCH THE TEST NUMBER DIRECTORY FILE (TDIR) FOR A 5PECI 

345 C COMBINATION OF TEST NUMBER, ENG INE NUMBER, OX I D I ZER PUMP NUMBER, FUE 

346 C PUMP NUMBER AND TEST DATE 

347 0 


349 

350 


U ubHbc 

C CALL STDIR (TEST, DATE, ENG, HPFP, HPOP, ISTAT) 

C 


D’ Line# 1 

351 C 

352 C 

353 C 

354 C 

355 C 

356 C 

357 C 

358 C 

359 C 

360 C 

361 C 

362 C 
b6b C 

364 

365 

•j>66 

367 

368 

369 

370 
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TEST - TEST NUMBER TO, BE SEARCHED FOR (6* CHACRACTERS) 

DATE - TEST DATE TO BE SEARCHED FOR (MO/DV/YR) 

ENG - ENGINE NUMBER TO BE SEARCHED FOR (10 CHARACTERS MAX. LEFT 
JUSTIFIED) 

HPFP - FUEL PUMP NUMBER TO BE SEARCHED FOR (10 CHARACTERS MAX. LEF 
JUSTIFIED) 

HPOP - OXIDIZER PUMP NUMBER TO BE SEARCHED FOR (10 CHARACTERS MAX. 
LEFT JUSTIFIED) 

ISTAT - SEARCH RESULTS 
1-SEARCH SUCCESSFUL 

O - NONE FOUND 
-1 - FILE READ ERROR 

SUBROUTINE STDI R (TEST, DATE, ENG, HPFP, HPOP, ISTAT) 

CHARACTER TEST*6, DATE*8, ENG*10, HPFP*10, HP0P*10 
CHARACTER TEST1*6, DATE1 *8, ENG1*10, HPFP1*10, HPOP 1*10 
CHARACTER TEST0*6, DATE0*8, ENG0*10, HPFP0*10, HP0P0*10, FD*4 
DATA TESTO, DATEO, ENGO, HPFPO, HPOPO /5*’ ’/ 

DATA FD/’ TDIR’ / 

IF ( ISTAT. NE. O) GO TO 1 


371 


LV=10 


372 


OPEN (LV, FILE=FD, IOSTAT=ISTAT) 

373 


IF ( ISTAT. NE.O) GC 

1 TO 999 

374- 

1 

CONTINUE 


375 


TEST 1 =’ O’ 


376 

5 

DATE1 = * O’ 


377 


ENGl^ 1 O’ 


378 


HPFP 1 =’ O’ 


379 

6 

HPQP1 =’ O’ 


380 

10 

READ ( 10, 1 1 , END=200 

, ERR=300) TEST 1 , DATE1 

381 


1ENG1 , HPFP 1 , HPOP 1 


382 

1 1 

FORMAT ( A6, IX, A8, 2X 

, 3(A10, 2X> ) 

383 

20 

CONTINUE 


384 


IF (TEST. NE. TESTO) 

GO TO 30 

385 

25 

CONTINUE 


386 


GO TO 40 


387 

30 

CONTINUE 


388 


I F ( TEST. NE. TEST 1 ) 

GO TO 10 

389 

40 

CONTINUE 


390 


IF (DATE. NE. DATEO) 

GO TO 50 

391 


GO TO 60 


392 

50 

CONTINUE 


393 


IF (DATE. NE. DATED 

GO TO 10 

394 

60 

CONTINUE 


395 


IF (ENG. NE. ENGO) GO 

TO 70 

396 


GO TO 80 


397 

70 

CONTINUE 


398 


IF (ENG. NE. ENG1 ) GO 

TO 10 

399 

80 

CONTINUE 


400 


IF (HPFP. NE. HPFPO) 

GO TO 90 

401 


GO TO 100 


402 

90 

CONTINUE 


403 


IF (mPFP. NE. HPFP 1 ) 

GO TO 10 

404 

1 00 

CONTINUE 


405 


IF (HPOP. NE. HPOPO) 

GO TO 110 

406 


GO TO 120 



* 


/ 

1 iU 

LUIM 1 1 IMUt 

408 


IF (HPOP. NE. HP0P1 ) GO TO 10 

409 

120 

1ST AT= 1 


D L l r»e# 
410 

1 

7 

411 


TEST=TEST 1 

412 

130 

DATE=DATE1 

413 


ENG=ENG1 

414 


hpfp=hpfpi 

415 


hpop=hpopi 

416 


RETURN 

417 

200 

CLOSE (10, IQSTAT=IERR> 

418 


ISTAT=0 

419 


RETURN 

420 

300 

ISTAT=-1 

421 


RETURN 

422 

999 

WRITE <*, 1000) ISTAT 

423 

1000 

FORMAT <’ FILE OPEN ERROR 

424 


STOP 

425 


END 
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14) 


Name 

Type 

Offset 

P 

DATE 

CHAR*8 

4 

* 

DATEO 

CHARTS 

14214 


DATE1 

CHAR*© 

14266 


ENG 

CHAR*1 0 

© 

* 

ENGO 

CHAR*1 0 

14222 


ENG1 

CHAR* 10 

14274 


FD 

CHAR*4 

14252 


HPFP 

CHAR* 10 

12 

* 

HPFPO 

CHAR* 10 

14232 


HPFP 1 

CHAR*10 

14284 


HP OP 

CHAR*1 0 

16 

* 

HPOPO 

CHAR*10 

14242 


HPOP 1 

CHAR* 1 0 

14294 


IERR 

INTEGER*4 

14328 


ISTAT 

INTEGER*4 

20 

* 

LV 

INTEGER*4 

14256 


TEST 

CHAR*6 

0 

* 

TESTO 

CHAR*© 

14208 


TEST 1 

CHAR*6 

14260 



Class 


426 C SEARCH THE TIME HISTORY FILES TO FIND THE MAGNITUDE AT A GIVEN TIME 

427 C USAGE: CALL STIMEH ( TEST, CH, T IME, MAG, ISTAT) 

428 C TEST-TEST NUMBER 

429 C CH-CHANNEL NUMBER 

430 C TIME-TIME FOR WHICH MAGNITUDE IS TO BE FOUND 

431 C MAG-MAGNITUDE AT TIME 

432 C ISTAT-CALL STATUS 

433 C — 3=T I ME EXCEEDS MAXIMUM TIME IN DATA FILE. MAXIMUM TIME AND 

434 C MAGNITUDE AT MAXIMUM IME ARE RETURNED IN TIME AND MAGNITUDE. 

435 C — 2=T I ME L.ESS THAN 0 

436 C — 1 =NORMAL RETURN 

437 C 1 —225— OPEN ERROR CODES 

438 C 

439 C NOT VERIFIED ON IBM 12/85 

440 C 

*+4 i 


SUBROUTINE STIMEH (TEST, CH, TIME, MAG, ISTAT) 


*+*+c. unhi^Mu i c.r\ i co i *ro, L^n ^o 

443 CHARACTER FD*1£, IFLT*1 

444 REAL MAG, MAGO, MAGI 
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D 


.ne# 

1 

445 

C 

446 


447 


448 

C 

449 

c 

450 


451 

10 

45£ 


453 


454 


455 


456 


457 


458 

£0 

459 


460 


461 


46£ 

50 

463 


464 


465 


466 


467 

70 

468 


469 

800 

470 


471 


473 


473 

850 

474 

875 

475 


476 


477 

1000 

478 


479 



09 
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**APPEARS TO BE MISSING FILTER DESCR 
WRITE (*,*) 5 ENTER IFLT, 1 CHAR, SUB STIMEH’ 

READ ( *, ’ (A1)M IFLT 
**ADDED l£/85 

WRITE (FD, 10) IFLT, TEST, CH 
FORMAT ( 1 T 1 , Al, A6, M » ,A3) 

LV= 1 £ 

I F < I STAT, NE. 0) RETURN 
0PEN(LV, FILE=FD, I 0STAT=I STAT) 

IF ( ISTAT. GT. 0) RETURN 
I ST AT =— 1 

READ (12, £0, END=800, ERR=850> XTIMEO, XMAGO 
FORMAT (£G13. 0) 

IF (XTIMEO. LT. TIME) GO TO 50 
I STAT =— £ 

GO TO 1000 

READ ( 1£, £0, END=SOO, ERR=850> XTIMEi, XMAG1 
IF ( XT I ME 1 . GE. TIME) GO TO 70 
X T I MEO= X T I ME 1 
XMAG0=XMAG1 
GO TO 50 

MAG=XMAGO+ (TIME- XTIMEO) * ( X MAGI— XMAGO) / ( XTIMEI -XTIMEO) 

GO TO 1000 
MAG—XMAG1 
TIME=XT I ME1 
I ST AT =— 3 
GO TO 1000 
WRITE (*, 875) FD 

FORMAT ( ’ READ ERROR ON FILE’,A1£) 

CLOSE (1£, IOSTAT==I ERR) 

STOP 

CLOSE (1£, I QSTAT=IERR) 

RETURN 

END 


3£ s 46 
0£/84 


Name 

Type 

CH 

CHAR*3 

FD 

CHAR*12 

I ERR 

INTEGER*^ 

IFLT 

CHAR*1 

ISTAT 

INTEGER*^ 

LV 

INTEGER*^ 

MAG- 

REAL 

MA GO 

REAL 

MAGI 

REAL 

TEST 

CHAR*6 

TI'vtE 

REAL 

XMAGO 

REAL 

XMAS 1 

REAL 

XTIMEO 

REAL 

XTIMEI 

REAL 


Offset P Class 

4 * 

14363 
1 445 £ 

1 436 £ 

16 * 

14396 
1£ * 

***** 

***** 

O * 

a * 

14404 
1 44EO 
14400 
14416 


480 CSPROG RDTH 


L, Ut)HtDC. S LHLL rtJU \ rt k I \ , Ufl, ITU, I il T IC., l T IHb, L-C.INJ, i'olHI ) 

482 C TEST-TEST NUMBER 

483 C CH— CHANNEL NUMBER 


D Line# 

484 

485 
46S 

487 

488 

489 

490 

491 

492 

493 

494 

495 

496 

497 

498 

499 

500 

501 

502 

503 

504 

505 

506 

507 

508 
• 509 

510 

511 

512 

513 

514 

515 

516 

517 
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C I FLT - FILTER CODE 
C TIME - TIME VECTOR 
C MAG-MAGNITUDE VECTOR 

C LEN - DIMENSION F TIME AND MAG ON CALL NUMBER OF DATA POINTS ON RETURN 
C ISTAT-CALL STATUS 
C 0=N0RMAL RETURN 

C — 1=NUMBER OF DATA POINTS GREATER THEN LEN 
C 1-225=0PEN ERROR CODES 

SUBROUTINE RDTH (TEST, CH, IFLT, TIME, MAG, LEN, ISTAT) 

CHARACTER TEST*6, CH*3 

CHARACTER FD*12 

REAL MAG (512) , T I ME (512) 

WRITE (FD, 10) IFLT, TEST, CH 
10 FORMAT ( 1 T’ A1 , A6, ’ . ’ , A3) 

LV=12 

IF (ISTAT. NE. O) RETURN 

OPEN (LV, FILE=FD, IOSTAT=ISTAT> 

IF ( ISTAT. NE. 0) RETURN 

LENM=LEN 

LEN=0 

15 LEN=LEN+1 

READ ( 12, 20, END=1000, ERR=850) TIME (LEN) , MAG (LEN) 

20 FORMAT (2G13. O) 

IF (LEN. LT. LEKlM) GO TO 15 ■ 

I STAT =— 1 
GO TO 1000 
850 WRITE ( *, 875) FD 

875 FORMAT ( ’ READ ERROR ON FILE’ , A12) 

CLOSE (12, IOSTAT=IERR) 

STOP 

1000 CLOSE ( 12, IOSTAT=IERR) 

LEN=LEN-1 

RETURN 

END 


Name 

T ype 

Offset 

P 

CH 

CHAR*3 

4 

* 

FD 

CHAR*1£ 

14456 


IERR 

INTEGER*^ 

14530 


IFLT 

INTEGER*^ 

8 

* 

ISTAT 

INTEGER*^ 

24 

* 

LEN 

INTEGERS 

20 

* 

LENM 

INTEGERS 

14490 


LV 

INTEGER*4 

14486 


MAG 

REAL 

16 

#* 

TEST 

CHAR*6 

0 

* 

TIME 

REAL 

12 

* 


Class 


518 CSPROG ERROR 

519 SUBROUTINE ERROR (M2SG, ISTAT) 

520 CHARACTER MES6*20 

521 WRITE(*, *> ’ SUBROUTINE ERROR’ 

522 WRITE <*, 1 ) MESG, ISTAT 

523 1 FORMAT ( ’ ERROR ’AEG,’ ISTAT= * 14) 


b I UP 

END 
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Narne Type Offset P Class 

I ST AT INTEGERS 4 * 

MESG CHAR*20 0 * 


526 C 

527 C$PR0G INTERP 

528 C 

529 C INTERPOLATION ROUTINE 

530 C 

531 SUBROUTINE INTERP < X, Y, XX, YY, SLOPE, LEN) 

532 DIMENSION X (500) , Y (500) 

533 IF (XX. GE. X ( 1 ) > GO TO ID 

534 YY=Y ( 1 ) 

535 RETURN 

536 10 XLST=X ( 1 ) 

537 YLST=Y ( 1 ) 

538 DO 20 J=2, LEN 

1 539 IF(X(J).GT. XX) GO TO 30 

1 540 XLST=X(J) 

1 541 20 YLST=Y(J> 

542 YY=Y (LEN) 

543 RETURN 

544 30 SLOPE=(Y(J)-YLST) / (X (J)-XLST) 

545 YY=YLST+SLOPE*(XX-XLST> 

546 RETURN 

547 END 


Name 

Type 

Offset 

P 

J 

INTEGER*4 

14570 


LEN 

INTEGER*4 

£0 

* 

SLOPE 

REAL 

16 

*■ 

X 

REAL 

0 

* 

XLST 

REAL 

1 4562 


XX 

REAL 

8 

* 

Y 

REAL 

4 

* 

YLST 

REAL 

14566 


YY 

REAL 

12 

*• 


548 CSPROG RDPSD 

549 C PURPOSE: READ THE CONTENTS OF A PSD DATTA FILE AND PUT 

550 C FREQ. IN X AND PUT MAG. IN Y AND PUT COMPOSITE 

551 C IN COMP 

552 C 


553 C 

554 C 

555 C 

556 C 

557 C 

558 C 

559 C 

560 C 


TEST - TEST NUMBER (1 WORD 6 CHARACTERS) 

CH - CHANNEL NO. 

I SLICE - SLICE CODE 

LEN - DIMENSION OF X AND Y WHEN CALLED - NUMBER OF DATA 
POINT ON RETURN 

' ! Y MUST BE DIMENSIONED ATLEAST LEN=10 IN MAIN PROGRAM 
Y (LEN J = LEN 

T ( LEN+2 ) = COMPOSITE ON PLOTS 


i 


OO X L_ 

5 62 c 

563 C 


1 

Y (LEN+4) 
Y (LEN+5) 


LUI'IHUDI it LUl'iHU I tU r fXUt»l UiDi 1 1 4LU 17M i M 

COMPOSITE COMPUTED FROM CONVERTED DATA 
NO. OF POINT NOT ASSIGNED DURING CONVERSION 
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564 C Y (LEN+6) = TIME SLICE 

565 C Y(LEN+7) = MAX . FREQUENCY 

566 C Y ( LEN+8) - Y(LEN+10) = SPARE 

567 C ISTAT - CALL STATUS 

568 C >0 OPEN ERROR CODES 

569 C , —500 C (—500— CLOSE ERROR CODES) 

570 C -1 READ ERROR 

571 C — £ DATA EXCEEDED DIMENSION 

57£ C LEN CHECK FAILED 

573 SUBROUTINE RDPSD (TEST, CH, ISLICE, X, Y, COMB, LEN, ISTAT) 

574 CHARACTER CH*3, TEST*6 

575 CHARACTER FD1£*12 

576 DIMENSION X < 500 ) , Y ( 500 ) 

577 WRITE (FD1£, 10) ISLICE, TEST, CH 

578 10 FORMAT <A£, A6, 1 . ’ A3) 

579 OPEN ( 12, FILE=FD1£, IOSTAT= I STAT) 

580 WRITE (5, ’ (A12) * )FD12 

581 WRITE ( *, *)FD1£ 

582 IF ( ISTAT- NE. 0) RETURN 

583 LENM=LEN+ 10 

584 LEN=0 

585 £0 LEN=LEN-+1 

586 READ <12, 30, END=50, ERR=40 ) Y ( LEN) 

587 IF (LEN. LT. LENM) GO TO £0 

588 30 FORMAT (G13. 6) 

589 I ST AT=— £ 

590 GO TO 60 

591 40 ISTAT=-1 

592 GO TO 60 

593 50 LENGLEN— 11 

594 COMB=Y (LEN+2) 

595 STEP=Y (LEN+7) /LEN 

596 DO 55 J= 1 , LEN 

1 597 55 X(J)=J*STEP 

598 C ?? IF (LEN. NE. I FI X (Y (LEN+1 ) ) ) ISTAT=-3 

599 60 CLOSE < 12, IOSTAT=IERR> 

600 IF ( I ERR- NE. 0) I STAT=- 1 ERR-500 

601 RETURN 

602 END 


Name Type 

CH CHAR*3 

COMB REAL 

FD12 CHAR* 12 
IERR INTEGER*4 

ISLICE INTEGER*4 

ISTAT INTEGER*4 

J INTEGER*4 

LEN I NTEGER*4 

LENM INTEGER*4 

STEP REAL 

TEST CHAR*6 
X REAL 

Y REAL 


Offset P Class 
4 * 


£0 

* 

14578 


14628 


a 

* 

28 

* 

1 4620 


24 

*■ 

14604 


1 46 1 6 


0 

* 

12 


16 

* 


SEARCH THE TEMPORARY PSD DIRECTORY FOR A GIVEN TEST NO 


SO 3 CPROG SPSD1 
£04 C PURPOSE: 


D Line# 

1 

605 

C 

606 

C 

607 

c 

608 


609 


610 


61 1 


612 


6 1 3 


614 


615 

820 

616 


617 

830 

618 


613 


620 


621 


622 

860 

623 


624 


625 


626 

1000 

627 

1 1 10 

628 


629 

1180 

630 

1 190 

631 


632 

2000 

633 

2010 

634 


635 
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AND CHENNEL NO. — RETURN WITH TIME SLICE 
MAX FREQ. AND SLICE CODE 

SUBROUTINE SPSD1 (TEST, CH, SLICE, FREQM, ISLICE, ISTAT, THRUST) 

CHARACTER TEST*6, TEST 1 *6, I SL I CE*2, CHX*3 
CHARACTER FD13*8, CH*3 
DATA FD13 /’ DPSD. TMP’ / 

IF ( ISTAT. EQ.O) GOTO 820 

OPEN < 13, FILE=FD13, I OSTAT= I STAT) 

IF ( ISTAT. NE. 0) GO TO 1180 

READ ( 13, 830, END=8S0, ERR=1000> TEST 1 , ISLICE, CHX, SLICE, 

1 FREQM, THRUST 

FORMAT (A6, IX, A2, 1 X, A3, 3 ( 1 X, G 13. 6 ) ) 

IF (TEST1. GT. TEST) GOTO 860 
IF (TEST1 . NE. TEST) GO TO 820 
IF (CHX. NE. CH) GO TO 820 
RETURN 
ISTAT—— 1 

CLOSE (13, IOSTAT=IERR) 

IF < IERR. NE. O) GO TO 2000 
RETURN 

WRITE <*, 1 110) 

FORMAT < T READ ERROR IN SPSD1M 
STOP 

WRITE(*, 1190) ISTAT 

FORMAT < 7 OPEN ERROR NO. ’13,’ ON FILE DSPD. TMPM 
STOP 

WRITE (*, 2010) IERR 

FORMAT (’ CLOSE ERROR NO. ’ 13’ ON FILE DPSD. TMP 1 ) 

STOP 

END 


Name 

Type 

CH 

CHAR*3 

CHX 

CHAR*3 

FD13 

CHAR*8 

FREQM 

REAL 

IERR 

INTEGER*4 

ISLICE 

CHAR*2 

ISTAT 

I NTEGER*4 

SLICE 

REAL 

TEST 

CHAR*6 

TEST 1 

CHAR*6 

THRUST 

REAL 


Offset P Class 

4 * 

14646 
14632 
12 * 

14678 
16 * 

20 * 

Q * 

O * 

14640 
24 * 


636 C 

637 C 

638 C 

639 C 

640 

641 

642 
£ 4^2 
644 


CREATE TEMP FILE FOR SELECTED CHANNEL NUMBER 
SEARCH DPSD FOR ALL ENTRIES OF MATCHING CHANNEL 
WRITE THEM TO DPSD. TMP 
SUBROUTINE DPSDTEMP ( TEST, CH, ISTAT) 

CHARACTER TE37*6, TEST 1 *6, ISLICE*2, CHX*3 
CHARACTER FD13*4, CH*3 
DATA PD13 /’ D^SD’ / 

GPEN ( 13, F I L£=FD 13, IOSTAT= ISTA T) 


to-H-D il" \ lb I H I « NC. m J DU IU ilOU 

646 I CHAN=0 

647 820 READ ( 1 3, 830, END=860, ERR=1000) TEST 1 , ISLICE, CHX, SLICE, 
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648 

649 

650 

651 

652 

653 

654 

655 

656 

657 

658 

659 

660 
661 
662 

663 

664 

665 

666 
667 
66 8 

669 

670 

671 

672 


1FREQM, THRUST 

830 FORMAT (A6, IX, A2, IX, A3, 3< IX, G13. 6) ) 

IF (CHX. NE. CH) GO TO 820 

IF (ICHAN. EG. 0) 0PEN(6, FILE = ’ DPSD. TMP’ > 

I CHAN= 1 

WRITE ( *, *) TEST 1 , CHX, ’ ’ , THRUST 

WRITE (6, 830) TEST 1 , ISLICE, CHX, SLICE, FREQM, THRUST 
GOTO 820 
860 CONTINUE 

IF ( ICHAN. EQ.O) WR I TE (*,*)’ NO PSD FILES FOR SELECTED CHAN’ 
CLOSE (13, IOSTAT=IERR ) 

IF ( IERR. NE. O) GO TO £000 
CLOSE (6, I QSTAT= IERR ) 

IF (IERR. NE. O) GOTO 2000 
RETURN 

1000 WRITE ( *, 1110) 

1110 FORMAT ( 1 READ ERROR IN DPSD’ ) 

STOP 

1180 WRITE (*, 1190) ISTAT 

1190 FORMAT ( 1 OPEN ERROR NO. ’13,’ ON FILE DSPD’ ) 

STOP 

2000 WRITE (*, 2010) IERR 

2010 FORMAT < * CLOSE ERROR NO. M3’ ON FILE DPSD OR . TMP’ ) 

STOP 

END 


Name 

Type 

CH 

CHAR*3 

CHX 

CHAR*3 

FD13 

CHAR*4 

FREQM 

REAL 

ICHAN 

INTEGER*4 

IERR 

INTEGER*4 

ISLICE 

CHARTS 

ISTAT 

INTEGER*4 

SLICE 

REAL 

TEST 

CHAR*6 

TEST 1 

CHAR*6 

THRUST 

REAL 


Offset P Class 

4 * 

14808 

14792 

14816 

14796 

14852 

14806 

8 * 

14812 

O * 

14800 
1 4820 


673 

674 

Name Type Size Class 


DPSDTE 

DSORT 

ERROR 

INTERP 

MAIN 

PEAK 

RDP SD 

RDTH 

SCHN 


SUBROUTINE 

SUBROUTINE 

SUBROUTINE 

SUBROUTINE 

PROGRAM 

SUBROUTINE 

SUBROUTINE 

SUBROUTINE 

SUBRCju T 1 |\E 


16 
;-86 
! s 46 
/84 


Of-OL/ 

3PSD1 
STD I R 


tsUtfKUU ■ i 'ic. 

SUBROUTINE 

SUBROUTINE 
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STIMEH SUBROUTINE 


Pass One No Errors Detected 
674 Source Lines 
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Section IV 


SUMMARY OF THE SSME DYNAMIC PRESSURE 
ENVIRONMENT AT RATED POWER LEVEL 

The appendices to this report include spectra for all the 
SSME tests, dynamic pressure measurements, and operating 
power levels provided by MSFC. A number of these 
measurements represent SSME operation over a limited power 
level range. However, the great majority of measurements 
include operation at RPL. In addition, an extensive data 
base of related vibration measurements has been developed at 
MSFC, representing SSME operation at RPL. 

Average spectra were therefore extracted from the data base, 
to summarize the dynamic pressure measurements obtained at 
RPL. The procedure applied is summarized as follows. 

1. Applying the data base sort routine, the spectra 
representing engine operation at RPL (99% - 101% RPL) were 
identified for each measurement. 

2. The spectra were plotted in a waterfall format to 
permit a qualitative evaluation of measurement trends and 
repeatabil ity . 

3. These spectrum families were reviewed in detail and, 
based on engineering judgement, clearly deviant spectra were 
deleted from subsequent analysis. (It should be noted that 
no spectra were deleted from the original data base.) 

4. Finally, the edited collection of spectra for each 
measurement was ensemble averaged using the statistical 
analysis subroutine, to obtain a single spectrum 


4-1 


representing each measurement at RPL. These results are 
presented as figures 4-1 through 4-27. 

The above described spectra should provide a valuable 
reference to summarize the SSME dynamic pressure environment 
under nominal operating conditions. In applying these plots 
as a "benchmark” , it should be noted that the spectra 
represent measurements with widely varying available sample 
size. For example, the auxilliary lox inlet pressure 
spectrum, figure 4-23, represents a single observation 
(PSD) , whereas the main combustion chamber hot gas inlet 
pressure spectrum, figure 4-8, is the average of 44 PSD's 
available in the data base at RPL. Therefore, considerably 
more variation is to be expected about the former 'average' 
spectrum than the latter. An interpolation scheme to 
estimate the pressure spectrum over a wide range of engine 
power levels is discussed in the next section. 
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Sample Size = 43 



FIGURE 4-5. SPECTRUM OF HPOP DS PR AT RPL 
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SECTION V 

SPECTRUM NORMALIZATION AND INTERPOLATION 

In the introduction to this report it was noted that the 
study was directed to accomplishment of three major 
objectives; data base generation, management software 
development, and development of techniques to estimate 
measurement spectral trends with engine operating 

conditions. This section summarizes approaches taken to 

accomplish the last objective. 

5.1 Dimensional Analysis 

At the outset of the study, it was anticipated that dynamic 
pressure measurement scaling would be accomplished through 
dimensional analysis based on governing fluid-dynamic and 
geometric parameters. A literature search was therefore 
initiated to define non-dimensional spectral relations for 
fluid/ structural systems similar to SSME components. Two 
such non-dimensional spectral models are as follows: 

Flow Induced Turbulence (1) (Duct Flow) 

S p = Q p (s)/P 2 v 3 d 

where Sp is the dimensionless pressure spectral density, and 
Qp = measured PSD 
P = local fluid density 
v = mean flow velocity 
d = characteristic length 

s= dimensionless frequency = fd/v (Strouhal number) 
f = frequency 


(1) Clinch, J.M. Prediction and Measurement of Vibrations 
Induced in Thin Walled Pipes by Internal Turbulent 
Waterflow. J. Sound Vib. , V. 12, No. 4, 1970. 
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figure 5-1 . Frequency spectrum of turbulent wall pres- 
sure field. (O), 267 in/sec; (A) 450 in/sec; (Q) 520 
in/sec 
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figure 5-2. Power spectral density of the high pressure 
oxidizer pump discharge pressure 5 seconds after SSME 
startup 



FIGURE 5-3- Nondimensionalized frequency spectrum for 
high pressure oxidizer pump discharge pressure, 5 
seconds after SSME startup 
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Blade Wake Excitation (2) (Pump Discharge) 


g(s, ) = Q p /Pc 2 D 3 M t 3 

where 

= Flow coefficient 
= Blade tip mach number 
s = Dimensionless frequency, f/n 

D = Blade Diameter 

n = Shaft speed 

Figure 5-1, from reference (1) , illustrates the non- 
dimensional wall pressure spectrum measured on a circular 
duct with internal water flow over a range of flow 
velocities. Collapse of the data is seen to be quite good. 
Figure 5-2 is the ordinary PSD of a dynamic pressure 

measurement. Dimensional normalization of this spectrum by 
the above technique is shown in figure 5-3 . The predominant 
spectral peaks appear related to turbopump rotational 
speeds. A normalization scheme was therefore developed to 
represent the spectra as a linear superposition of simple 
non-dimensional sources. 

As sufficient data base measurements became available a 

detailed evaluation of spectral trends with engine power 
level was performed. (It should be noted that the non- 
dimens ionali zing parameters such as flow velocity and 

turbopump blade tip speed vary in proportion to engine power 
level) . Since well defined spectral trends were not 

indicated in general, it became clear that the measured flow 
dynamic interactions are too complex to be adequately 
represented by a superposition of simple sources. 


(2) Maling, G. C. Jr. Dimensional Analysis of Blower Niose. 
J. Acoust. Soc. Am., V. 35 No. 10, 1963. 
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Dimensional analysis was therefore abandoned in favor of a 
more general two dimensional interpolation procedure. 
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5.2 Interpolation Bv Singular Value Decomposition 


The method of Singular Value Decomposition is a mathematical 
tool for the decomposition of general rectangular matrices. 
By this technique, a family of curves, such as the power 
spectra of dynamic pressure measured at a number of power 
levels can be reduced to a superposition of basic component 
patterns (eigenvectors) and associated weighting 
coefficients. The predominant normalized component patterns 
yield a good approximation of the given spectra. Thus, the 
difficult problem of interpolating (two dimensional) over 
the observed spectra can be reduced to one dimensional 
polynominal interpolation of the derived weighting 
coefficients. 

It may be instructive to observe that the SVD technique is 
mathematically akin to the normal mode method in structural 
dynamics. In the latter application, the complex deflection 
patterns of a structure are expressed as a superposition of 
normal modes, or eigenvectors. The associated weighting 
coefficients (eigenvalues) are the system natural 
frequencies. (An excellent review of the SVD technique, 
upon which this discussion and the resulting software is 
based, is given by S. Tominaga in the IEEE Transactions on 
Acoustics, Speech and Signal Processing, Vol. ASSP-29, No. 
3, June, 1981.) The technique is best illustrated by 
example. 

Figure 5-4 illustrates the ensemble average spectra of the 
high pressure oxidizer pump discharge pressure observed at 
11 engine power levels. Based on this set of data we wish 
to extract normalized spectrum functions and associated 
coefficients to represent the curves over a continuous range 
of power levels. Each PSD may be represented as a vector 
p, (n x 1) , representing the spectral values at n 
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frequencies. The group of spectra may thus be written as a 
matrix P, (n x m) , representing the m spectra at n 
frequencies. (In the case at hand this is a 400 x 11 
matrix.) If the rank of P is m, it may be expressed as the 
sum of m matrices of rank one, 

P = r l n l e T l + . . . + r m n* e T m 

where n^, e^ are the orthonormal vectors of PP^ and P^P, 
respectively, and the singular values r^ are the square 
roots of the eigenvalues of PP T . 

Now, each of the eleven spectra may be uniquely defined as a 
linear combination of the fundamental component waves n^, in 
the form 

Pi = c il n l + •••• + c im n m 

where c^j = and ej is defined to be [e-^, ...e m j] T . 

The SVD technique provides a method of synthesizing a set of 
observed curves. For purposes of prediction it is then only 
necessary to interpolate over the coefficients c^j , as a 
function of engine power level, for each fundamental 
component pattern. 

The software was developed for the implementation of the SVD 
technique and interpolation scheme. A program listing of 
the routines is given in Table 5-1. Figure 5-5 illustrates 
the fundamental component patterns associated with the three 
largest singular values for the spectra in Figure 5-4. 

The above analysis indicates that all component patterns be 
utilized to synthesize the family of curves. This provides 
a unique representation. In practice, however, convergence 
is quite rapid. This is again similar to the normal mode 
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FIGURE 5-k. ENSEMBLE AVERAGE SPECTRA AS A FUNCTION OF ENGINE POWER 
LEVEL (HPOP DS PR) 


NORMALIZED AMPLITUDE 




0 5000 10,000 
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FIGURE 5 - 5 . NORMALIZED SPECTRUM FUNCTIONS (EIGENVECTORS) 

5-9 


* w 


■J-P 

<T\ 

O 


‘V*. 


L .! . 

1 feV. J .1 


■f. ft 

-<J£j 1 ■!%;=■ 




TBS J i tr 


\ 


»- ■. i »- * >. ‘V. 

5 i i i ^ u sVTS^.c^ ^ 

< v! 1 ? -ai.S* 


Q- 

<*e 

o 

vO 



\ 

l 

i- 

i 

! 

1 

i 


I— 

i 




CSS 


a: 

< 


co 

>- 
00 

0 
2 
LU 

oc: lu 

I- 2 

LU 

—i a: 
< o 
a: 2 

h— ‘ — 
CJ 

LU U 

01 
CO 

o 

LU 


t 

1 

f— 

i 

3 

i 


CO 

t: 


o. 


< 5 * 

ie\ 


s ^2 

cc 

rCTT W 

» t- 


NO 

I 

LT\ 


ce: 

ZD 

o 


C^l 


to 

o 

C- 

n 

o 

CD 


y 

_ j 
< 
> 


( ZH/ ISJ ) (ISi DO! 

I 




5-10 


method, wherein the ’deflection of a beam can generally be 
adequately represented by the superposition of a limited 
number of modes. 

For the example under discussion, the eleven computed 
singular values (we have eleven PSD's) in descending order 
are: 105.3, 9.8, 6.5, 4.4, 3.9, 3.4, 3.3, 2.9, 2.7, 2.3 and 
.00002. The relative contribution of each fundamental 
pattern may be shown to be proportional to its associated 
singular value. It is therefore clear that a superposition 
of the first three fundamental patterns should provide a 
good approximation to the given family of spectra. This was 
done for the cited example, and an interpolated library of 
pressure spectra was computed at 2.5% power levels over the 
range of observed test values. These results are shown in 
Figure 5-6, and indicate the spectral trend of dynamic 
pressure with engine power level. 

5.3 An Interpolated Library of Spectra 
for Dynamic Pressure Measurements 

To indicate the trend in amplitude/ frequency of the SSME 
dynamic pressure measurements an interpolated library of 
spectra was derived, by the method just discussed. The 
analysis procedure was performed as follows: 

1. The data base spectra for each measurement 
location were plotted, and any obvious outliers 
discarded from subsequent analysis. 

2 . The remaining data sets were segregated by power 
level. The range and number of power levels available 
varied by measurement location, based on the specific 
measurements and test conditions provided by MSFC. 
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3 . The spectra for each measurement/power level were 
ensemble averaged to obtain a family of spectra as 
shown in the example. Figure 5-4. 

4. An SVD analysis was performed on each set of 
curves to extract intrinsic component patterns and 
associated singular values. The number of derived 
patterns/coefficients varied according to the number of 
test power levels available for each measurement. 

5. Finally, interpolation was applied to the above 
derived coefficients, and plots generated representing 
spectral variation in 2.5% power level increments over 
the available data range. 

The results of this exercise are illustrated in the 
remaining figures of this report. These figures provide a 
qualitative indication of spectrum variation with engine 
power level. A quantitative estimate of spectral amplitudes 
and frequencies may be obtained by cross reference to the 
average spectra representing operation at RPL, as given in 
Section IV of this report. 

It will be noted that several of the data base measurements 
are not included in the above discussed spectrum library. 
This is essentially due to lack of sufficient data; either a 
very limited number of spectra and/or a very limited test 
power level range. In either case, one is referred to the 
average spectra presented in Section IV or the data base 
plots provided in the Appendices for spectral definition of 
these measurements. 
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5.4 A Computer Routine for Estimating 
Dynamic Pressure Environments 


The interpolated spectrum libraries presented in the 
previous discussion provide a valuable indication of 
spectral trends with engine power level. However, as noted, 
the results are qualitative. A computer routine was 
therefore developed to generate an estimated spectrum for a 
user defined dynamic pressure measurement and engine power 
level . 

This program utilizes the eigenvectors and singular values 
obtained by the SVD technique as discussed in Section 5-3 
for each measurement, and is therefore totally indepedent of 
the complete data base. The routine was designed for the 
engineer responsible for dynamic evaluations and requires no 
knowledge of the linear algebra implied by SVD. 

A simple menu driven format is provided, wherein the user is 
querried as to what measurement description and power level 
a PSD estimate is desired for. (The user is also asked how 
many eigenvectors (fundamental SVD patterns) he/ she desires 
to use in the spectrum estimate. To answer this question 
the user has the option to review the singular value 
magnitudes associated with each pattern. However, through 
development of the SVD analysis and the pressure spectrum 
libraries, experience indicates that a "three mode" 
approximation should be satisfactory and efficient.) A 
program listing of the spectrum estimation routine is given 
in Table 5-1. The program is operational on any 
PC/compatible computer. 
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FIGURE 5-7. SPECTRAL TREND OF LPFP IN PR (60 - |00% RPL, 2.5 * INCREMENTS) 
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FIGURE 5-8. SPECTRAL TREND OF LPOP DS PR (60 - 100% RPL , 2.5% INCREMENTS) 
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FIGURE 5-9. SPECTRAL TREND HPOP BAL CAV PR2 (60 - 109$ RPL) 
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FIGURE 5-10. SPECTRAL TREND OF HPOP BAL CAM PR1 (60 - 109% RPL) 
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FIGURE 5-12. SPECTRAL TREND LOX IN DUCT PR2 (60 - 100% RPL) 
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FIGURE 5-U*. SPECTRAL TREND MCCHOT GAS IN PR (59 - 100% RPL) 
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FIGURE 5-16. SPECTRAL TREND LPFP DS PR (60 - 100% RPL, 2 .5% INCREMENTS) 
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FIGURE 5-17. SPECTRAL TREND PBP FUEL SUB PR (60 - 30 % RPL) 
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FIGURE 5-18. SPECTRAL TREND FPB FUEL MAN PR (60 - 100% RPL) 
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FIGURE 5-19. SPECTRAL TREND LPOT TURB DR PR (60 - IOOI RPL) 
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FIGURE 5-20. SPECTRAL TREND LOX HX OUT PR (60 - 90% RPL) 


100 % 



5-28 


FIGURE 5-21. SPECTRAL TREND HPFP BAL CAV PR (65 - 100% RPL) 
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FIGURE 5-22. SPECTRAL TREND MCC FUEL l NJ PR (65 - 109% RPL) 
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FIGURE 5-23. SPECTRAL TREND FPB PC (65 - 109% RPL) 
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FIGURE 5-24. SPECTRAL TREND MCC FUEL I N J PR DC (65 - 109^ RPL) 
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FIGURE 5-26. SPECTRAL TREND OPF PC PR DC (65 - 109% RPL) 


TABLE 5-1. A SPECTRUM INTERPOLATION PROGRAM 
BASED ON SINGULAR VALUE DECOMPOSITION 


Page 1 
01 -SB- 86 


Line# 

1 


4 

5 

6 
7 

a 

9 

10 
i i 
12 

13 

14 

15 

16 
17 
13 
19 
£0 
£1 


£4 

£5 

£6 

£7 

£8 

£9 

30 

31 
3£ 

33 

34 

35 

36 

37 
36 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 


09 : 36 : 45 

1 7 Microsoft F0RTRAN77 03.20 02/84 

C PROGRAM SVDNA5A 

C THIS PROGRAM GENERATE THE PSD OF SSME PRESSURE MEASUREMENTS 

C FROM 31 DIFFERENT LOCATIONS. THE ORIGINAL PSD HAVE BEEN 

C TRANSFORMED TO A FEW BASIC WAVEFORMS BY USING SINGULAR 

C VALUE DECOMPOSITION METHOD. THESE BASIC WAVE FORMS WERE 

C USED IN THIS PROGRAM TO ESTIMATE THE PSD AT ANY POWER LEVEL 

C WITHIN CERTAIN RANGE. 

REAl_*4 XXI ( 100) , XX2 ( 100) , WK1 ( 100) , WK2 (100) 

DIMENSION ID (34) , NN (34) 

REAL*4 CC ( 16, 16) , DD ( 16, 16) , U (400, 12) , VV (400, 12) , PSD (400) 
CHARACTER*3 FNAME ( 34 ) 

CHARACTER* 1 0 F1,F2,F3 
CHARACTERS 1 F11,F12,F13 
CHARACTER* 1 6 F4, F5 
WRITE (*, 666) 

666 FORMAT ( ’ THIS PROGRAM GENERATES THE PSD OF PRESSURE MEASURE’, 

1 ’ MENT AT ANY POWER LEVEL’,/,’ ( WITHIN CERTAIN RANGE ) BY ’ , 

1 ’USING SINGULAR VALUE DECOMPOSITION METHOD’) 

WRITE (*, ’ (A) ’ ) ’ OUTPUT PSD HAS 400 LINE SPECTRA’ 

OPEN ( 1, FILE"* SVDNASA. DAT’ , STATUS=’ OLD’.) 

DO 12 1=1,34 

READ ( 1 , 60 1 ) FNAME ( I ) , ID ( I ) , NN ( I ) 

12 CONTINUE 

601 F0RMAT(A3, IX, II, IX, 12) 

CLOSE (1) 

M=400 

123 CALL CHANNEL 

WRITE (*, ’ (A) ’ ) ’ ENTER CHANNEL #’ 

READ (*, *) NCH 

IF ( ID (NCH) . EQ. 9) WRITE (*, ’ (A) ’ ) ’ NO DATA 1 ’ 

IF (ID (NCH). EQ. 9) GOTO 123 

CCCCCCCCCCCCCCCCCCCCCCCCCCC I D=2 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
IF ( ID (NCH) . EQ. 2) THEN 

WRITE (*,’ (A)’)’ ENTER 1 FOR MAX FREQ = 20000 HZ ’ 

IF (NCH. NE. 7)WRITE(*, ’ (A)’ )’ 2 FOR MAX FREQ = 40000 HZ ’ 

IF (NCH. EQ. 7) WRITE (*, ’ (A) ’ ) ’ 2 FOR MAX FREQ = 5000 HZ’ 

READ ( *, *) IF12 

I F ( I F 1 2. EQ. 1 ) THEN 
WRITE (FI 1, 602) FNAME (NCH) 

WRITE (F12, 603) FNAME (NCH) 

WRITE (F13, 604) FNAME (NCH) 

ELSE 

WR I TE ( FI 1 , 605) FNAME (NCH) 

WRITE (F12, 606) FNAME (NCH) 

WR I TE ( F 1 3, 607) FNAME (NCH) 

END IF 

OPEN ( 1 , F I LE=F 11, STATUS=’ OLD’ , ACCESS=’ DIRECT’ , RECL=4*M) 

□ PEN (2, FILE=F12, STATUS=’ OLD’ ) 

OPEN (3, FILE=F13, STATUS=’ OLD’ ) 

WRITE (*, ’ (A) ’ ) Fi 1 , F12, F13 
END I ~ 



54 


56 

57 
53 
59 


CCCCCCCCCCCCCCCCCCCCCCCCCCC I D=0, 


ccccccccccccccccccccccccccccc 


IF ( ID (NCH) . LE. 1 ) THEN 

IF (NCH. NE. 31 ■ AND. NCH. NE. 30. AND. NCH. NE. 9) i Hc.N 
WRITE!*,’ (A)’)’ MAX FREQ. OF OUTPUT PSD = 20000. HZ’ 


’ (A) 1 ) ’ MAX r'KzO. Or UU i : J U ; Ph>D = 40000. HZ’ 

END I F 

WR I TE ( F 1 , 60S ) FNAME ( NCH ) 


5-34 


D 


1 

i 

i 

i 


i 

l 

i 


L i rie# 
60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 

32 

83 

84 

85 

86 

87 

88 

89 

90 

91 
98 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 
109 
1 10 
ill 
lie! 
1 t3 
t 1 4 
1 15 
1 IS 
1 17 
118 


Page <=: 
01-28-86 
09 : 36 : 45 

1 7 Microsoft F0RTRAN77 03.20 02/84 

WRITE (F2, 609) FNAME (NCH) 

WRITE (F3, 610) FNAME (NCH) 

□ PEN ( 1, FILE=F1, STATUS=’ OLD’ , ACCESS=’ DIRECT’ , RECL=4*M> 

IF ( ID <NCH) . EQ. 1 ) OPEN (2, F ILE=F2, STATUS=’ OLD’ ) 



OPEN (3, 

t- ILt=t-3, 

ST ATUS=’ 


ENDIF 



602 

FORMAT (’ 

PSD’ , A3, 

’ A. EVT’ ) 

603 

FORMAT (’ 

PSD’ , A3, 

’ A. EVl_’ > 

604 

FORMAT (’ 

PSD’ , A3, 

’ A. PRO’ ) 

605 

FORMAT (’ 

PSD’ , A3, 

’ B. EVT’ ) 

606 

FORMAT (’ 

PSD’ , A3, 

’ B. EVL’ ) 

607 

FORMAT (’ 

PSD’ , A3, 

’ B. PRO’ ) 

608 

FORMAT (’ 

PSD’ , A3, 

’ . EVT’ ) 

609 

FORMAT (’ 

PSD’ , A3, 

’ . EVL’ ) 

610 

FORMAT (’ 

PSD’ , A3, 

’ . PRO’ ) 


M=400 
N=NN (NCH) 

READ(3, *) (XXI (I ) , 1=2, N+l) 

CLOSE (3) 

DO 22 J=1,N 

READ ( 1 , REC=J) (U(I, J) , 1 = 1, M) 

22 CONTINUE 

DO 1 J=1,N 

IF ( ID (NCH) . NE. 0) READ (2, *) <CC(I,J), 1=1, N) 

1 CONTINUE 

CLOSE ( 1 ) 

IF (ID (NCH) . NE. 0) CLOSE (2> 

321 N0UT=1 

WRITE (*,’ (A)’)’ ENTER FILE NAME OF OUTPUT PSD’ 

READ ( *, ’ (A) ’ )F5. 

OPEN (5, FILE=F5, STATUS=’ NEW’ ) 

CCCCCCCCCCCCCCCCCCCCCCCC ID=0 CCCCCCCCCCCCCCCCCCCCCCCCCCC 
IF ( ID (NCH) . EQ. 0) THEN 

621 FORMAT ( ’ PSD ONLY AVAILABLE AT POWER LEVEL: ’ , 1 0 ( F7. 2, ’ •/. ’)) 
WRITE (*, 621 ) (XXI ( I ) , 1=2, N+l ) 

WRITE (*, ’ (A) ’ ) ’ ENTER DESIRED POWER (IN '/.) ’ 

READ (•*,*> PWL 
AMIN=200 
M I N=2 

DO 691 1=2, N+l 

TE=ABS (PWL-XX1 ( I ) ) y 

IF (TE. LE. AMIN) THEN ' 

AMIN=TE 
M I N= 1 - 1 
ENDIF 

691 CONTINUE 

WRITE (5,*) ( U ( I , MIN) , 1=1, M) 

CLOSE (5) 

GOTO 1986 
ENDIF 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
IF ( ID (NCH) . EQ. 1 ) THEN 
WRITE (*,611) XXI (2) , XXI (N+l ) 

611 FORMAT (’ AVAILABLE POWER LEVEL = ’ , FxO. 4, ’ % TO ’ , F10. 4, ’ '/.’ ) 
WRITE (*, ’ (A) ’ ) ’ ENTER DESIRED POWER LEVEL (%) ’ 

READ ( *, *) XX2 ( 1 ) 

ELSE 

WRITE (*, ’ (A) ’ ) ’ ENTER INPUT POWER LEVEL 
READ (*, ’ (A) ’ ) F 4 


F ILNAM’ 


1 


Pag 

01 - 


D Line# 
1 19 
120 
121 
122 
1 £3 

124 

125 

126 
127 
12Q 

129 

130 

1 131 

2 1 32 

2 133 

134 

1 135 

2 1 36 

1 137 

1 138 

1 139 

1 140 

1 141 

1 142 

1 143 

1 144 

2 145 

1 146 

147 

148 

1 149 

2 1 50 

2 151 

3 152 

1 jJi 

1 154 

1 155 

2 156 

1 157 

1 158 

1 159 

160 

161 

162 

163 

164 

165 

166 

167 

168 
169 


09 


61 1; 


1 1 


65 


64 

66 


7 Microsoft hGRjRAN77 V3. 20 

OPEN <4, FILE=F4, STATUS^ OLD’ ) 

READ (4, ■*) ( X X2 ( I ) , 1=1, NOUT) 

CLOSE (4) 

END IF 
IFLAR=2 
WRITE ( *, 61£)N 
READ <*, *) PIMM 

FORMAT <’ ENTER # OF BASIC WAVEFORMS TO BE USED <LE. ’ , 13, ’ ) * ) 
IFL=1 

XXI < 1 ) =XX 1 (2)— (XXI (3)— XXI (2) ) . 

XXI <N-*-2) =XXKN+1) + (XX1(N+1)-XX1(N)) 

DO 11 1=1, M 

DO 11 J=1,MMM 
W(I, J)=U<I, J) 

CONTINUE 
DO 66 1=1, MMM 
DO 65 J=1,N 
WK1 <J+1)=CC(I, J) 

WK1 ( 1 ) =WK 1 (2) + (WKl <2> — WK 1 (3) ) 

WK1 < N+2) =WK 1 (N+l ) 4- (WK1 <N+1 ) -WK1 (N) ) 

I F ( IFLAR. EQ. 1 ) THEN 

CALL LAR (XXI, WK1, N+2, NOUT, XX2, WK2) 

ELSE 

CALL LINE <XX1, WK1, N+2, NOUT, XX2, WK2 ) 

END IF 

DO 64 J= 1 , NOUT 
DO ( I , J)=WK2(J) 

CONTINUE 


C*****************************-**-* 

DO 102 11=1, NOUT 
DO 101 1=1, M 
PSD ( I ) =0. 

DO 101 K= 1 , MMM 

PSD ( I ) =PSD < I ) +W < I , K)*DD<K, II) 

101 CONTINUE 

IF(IFL. EQ, 1 ) THEN 
DO 104 LL= 1 , M 

104 PSD < LL ) =10. ** ( PSD ( LL ) ) 

ENDIF 

WRITE (5, 103) ( PSD < I ) , 1 = 1, M) 

102 CONTINUE 

103 FORMAT ( IX, 5E15. 7) 

CLOSE <5) 

1986 WRITE (*, ’ (A) ’ ) ’ ENTER 1 FOR SELECTING ANOTHER POWER LEVEL’ 
WRITE <*, 1 (A)’)’ 2 FOR SELECTING ANOTHER CHANNEL’ 

WRITE (*, ’ (A)’)’ * ELSE FOR STOP PROGROM’ 

READ ( *, *) I FCG 

I F ( I FCG. EQ. 2) GOTO 123 

IFCIFCG.EG. DGOTO 321 

STOP 

END 


02/84 


Name 

i yoe 

Otfset P Class 

ABS 


, INTRINSIC 

N 

^£4i_ 

hM-566 

OC 

REAL 

4:87m- 

DD 

*£AL 

42398 

F 1 

Ci-AR * 1 0 

44261 


Oi "J ro 


D Line# 1 7 


FI 1 

CHAR* 11 

44538 

FI 5 

CHAR*! 1 

44539 

r 13 

CHAR* 1 1 

44550 

F5 

CHAR* 10 

44571 

F3 

CHAR* 10 

44581 

F4 

CHAR* 1 6 

44638 

F5 

CHAR* 1 6 

44490 

FNAME 

CHAR*3 

43933 

I 

INTEGER*4 

44196 

ID 

INTEGER*4 

1603 

I FI 5 

INTEGER*4 

44354 

IFCG 

INTEGER*4 

44780 

IFL 

INTEG£R*4 

44716 

IFLAR 

INTEGER*4 

44655 

II 

I NTEGER*4 

44740 

J 

I NTEGER*4 

44474 

K 

INTEGER*4 

44755 

L.L 

INTEGER*4 

44760 

M 

INTEGER*4 

44316 

MIN 

INTEGER*4 

44570 

MMM 

II\1TEGER*4 

44656 

N 

INTEGER*4 

44466 

NCH 

l!MTEGER*4 

44350 

NN 

INTEGER*4 

1736 

NOUT 

INTEGER*4 

44486 

PSD 

REAL 

40374 

PWL 

REAL 

44563 

TE 

REAL 

44578 

U . 

REAL 

1874 

VV 

REAL 

31074 

WK1 

REAL 

803 

WK3 

REAL 

1303 

XXI 

REAL 

5 

xxs 

REAL 

403 
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1 

1 


i 


1 


1 


170 

171 
175 

173 

174 

175 

176 

177 

178 

179 

180 
181 
183 

183 

184 
165 
186 
187 
138 
169 
190 
19 l 

t 9 a 


c*****************-****************************-** 


SUBROUTINE LAR ( XX 1 , YY 1 , N IN, NOUT, XX3, YY3> 

REAL*4 XXI <1>, YY1 (1), X X5 < 1 ) , YY2 < 1 ) , CC < 50 ) , FF < 50) 

DO 3 1=1, NIN 

CC(I>=1. 00 

DO 4 J=l, NIN 

IF ( J. EQ. I) GOTO 4 

CC(I)=CC(I)*(X,X1 (I) -XXI (J) > 

CONTINUE 

CC ( I > =YY 1 <I)/CC<I) 

CONTINUE 

DO 5 1=1, NOUT 

A=XX5(I) 

DO 6 J=1,NIN 

FF ( J) =1 . 00 

DO 7 K=l, NIN 

I F ( K - EQ. J ) GOTO 7 

FF ( J ) -FF (J) * (A-XX1 <K) ) 

CONTINUE 

FF<J)=FF<J)*CC<J> 

CONTINUE 


VYU ( I ; 


t ) 


xJ U 


S '_=1 , N 1 N 


7 


D Line# 

‘ 1 

7 



2 193 


YY2 

< I ) =YY2 < I ) • 

2 194 

8 

CONTINUE 


1 195 

5 

CONTINUE 


196 

RETURN 


197 


END 



Name 

Type 


Offset 

P 

A 

REAL 


45204 


CC 

REAL 


44784 


FF 

REAL 


44984 


I 

INTEGERS 

45184 


J 

INTEGER*4 

45192 


K 

INTEGER*4 

45212 


L 

INTEGER*4 

45220 


NIN 

INTEGERS 

a 

* 

NOUT 

I NTEGER*4 

12 

* 

XXI 

REAL 


O 

*■ 

XX2 

REAL 


16 

-* 

YY1 

REAL 


4 

* 

YY2 

REAL 


20 

* 
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19Q 

199 

200 

1 201 
2 202 
2 203 

2 204 

2 205 

2 206 1 

1 207 93 

208 
209 


SUBROUTINE LI NE ( XX 1 , YY1 , NI N, NOUT, XX2, YY2 ) 

REAL*4 XXI < 1 ) , YY1 < 1) , XX2 < 1 > , YY2 ( 1 > 

DO 99 1=1, NOUT 
DO 1 J=1,NIN-1 

IF (XX2 ( I ) . GE. XXI < J) . AND- XX2 < I > „ LE. XX 1 < J+ 1 ) > THEN 

YY2 ( I) =SSS (XXI ( J) , XXI (J+l) , YY1 ( J) , YY1 ( J+l ) , X X2 < I ) ) 

GOTO 99 

END IF 

CONTINUE 

CONTINUE 

RETURN 

END 


Name Type 

I INTEGER*4 

J I NTEGER*4 

NIN INTEGERS 

NOUT INTEGER*4 

XXI REAL 

XX2 REAL 

YY1 REAL 

YY2 REAL 


Offset P Class 

45228 

45236 

a * 

12 * 
o * 

16 * 

4 * 

20 * 


210 
21 1 
212 
213 


FUNCTION SSS<X1, X2, Yl, Y2, XX) 
SSS=Y1+ (XX-Xl ) * <Y2-Y1 ) / (X2-X1 ) 
RETURN 
END 


Name 

Type 

Offset 

p 

XI 

REAL 

0 

* 

X2 

REAi, 

4 

*- 

XX 

REAL 

16 

* 

Yl 

REAL 

a 

* 


C 1 a S5 


1 


Page 6 
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Y2 REAL 12 * ' 


214 SUBROUTINE CHANNEL 

215 WRITE (*, ’ (A) ’ ) ’ 1. THRUST (NO DATA) 12. LPFP DS PR 

216 1 23. FPB PC’ 

217 WRITE <*, ’ < A) ’ ) ’ 2. LPFP IN PR 13. PBP FUEL SUP PR 

218 1 24. MCC FUEL INJ PR DC 1 

219 WRITE (*, ’ (A) ’ ) * 3. LPOP DS PR 14. FPB FUEL MAN PR 

220 1 25. OPB PL*’ 

221 WRITE (*, ’ (A) ’ ) ’ 4. HPOP BAL CAV PR2 15. LPOT TURB DR PR 

222 1 26. HPOP IN PR’ 

223 WRITE (*, ’ < A) ’ ) ’ 5. HPOP BAL CAV PR1 16. LOX HX OUT PR 

224 1 27. AUX LX I PR*’ 

225 WRITE <*, ’ ( A) ’ ) ’ 6. HPOP DS PR 17. OPB PC PRR* 

226 1 28. HI LX I PR*’ 

227 WRITE <*, ’ <A) ’ ) ’ 7. LOX IN DUCT PR2 18. HEX OUTLET* 

228 1 29. HOT GAS MAN PR2*’ 

229 WRITE (*, ’ (A) ’ ) ’ 8. HPOP DISC PR 19. HPFP DS PR* 

230 1 30. MCC LONG ACCEL 1*’ 

231 WRITE <*, ’ (A) ’ ) ’ 9. MCC HOT GAS IN PR 20. HOT GAS MAN PR3* 

232 1 31. MCC IN PR*’ 

233 WRITE <*, ’ < A) ’ ) ’ 10. PBP DS PR 21. HPFP BAL CAV PR 

234 1 34. OPB PC PR DC’ 

235 WRITE <*, ’ (A)’ ) ’ 11. CHAMBER PR* 22. MCC FUEL INJ PR ’ 

236 WRITE (*,666) 

237 WRITE <*, 667) 

238 666 FORMAT ( ’ NOTE: CHANNEL DESCRIPTOR ENDED WITH ,, * H MEANS’, 

239 1 ’ SVD DATA IS NOT AVAILABLE’) 

240 667 FORMAT ( ’ PSD ONLY AVAILABLE AT SOME FIXED POWER LEVEL’) 

241 RETURN 

242 END 


Name Type Offset P Class 


Name 

Type 

Size Class 

CHANNE 


SUBROUTINE 

LAR 


SUBROUTINE 

LINE 


SUBROUTINE 

MAIN 


PROGRAM 

SSS 

REAL 

FUNCTION 


Pass One 


No terrors Detected 
242 Source Lines 


D L 1 ne# 
1 
£ 

3 

4 


1 

P 


1 

i 

S'. 

i 

i 


c 

c 

c 


6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 
IS 

19 

20 
21 
22 

23 

24 

25 

26 
27 
£8 

29 

30 1 

31 

32 

33 

34 78 

35 79 

36 

37 

38 

39 

40 

41 

42 502 

43 

44 501 

45 

46 

47 

48 

49 504 

50 

5 1 503 

52 

53 

54 

55 777 

56 

57 

sa 

59 333 
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THIS PROGRAM USES SINGULATR VALUE DECOMPOSITION METHOD TO * 

GENERATE THE BASIC WAVEFORMS AND ITS SINGULAR VALUES IN ORDER 
TO ESTABISH THE WHOLE FAMILY OF PSD FROM THE GIVEN MEASURED 
PSD AT SOME CERTAIN PARAMETER ( e. g. POWER LEVEL ) 

DIMENSION IJK<12), WKC12) 

REAL*4 CC < 1 1 , 1 1 ) , DD < 1 1 , 1 1 ) 

REAL*4 A <M, N) , U <M, N> , V <M, N) , SIGMA <M) , WORK(M) 

INPUT # OF DATA IN PSD MUST EQUAL M 
INPUT # OF PSD COULD BE SMALL THEN N 

REAL*4 A <400, 11), U <400, 1 1 ) , V <400, 11), SIGMA <400) , WORK <400 
CHARACTER* 16 F1,F2,F3 

WRITE <*,’ <A) ’ ) ’ ENTER # OF DATA IN PSD < MUST EQUAL 400)’ 

READ <*, *) M 

WRITE <*, ’ < A) ’ ) ’ ENTER # OF PSD <LE. 11)’ 

READ <*, * ) N 
WRITE <*, * <A) 

READ <*, ’ 

WR I TE < *, 

READ <*, * 

WRITE < *, 

READ <*, ’ 

WRITE <*, 

READ <*, *) IFST 
WRITE <*, » <A> » ) ’ 

READ <*,*) LOGY 

OPEN < 1, FILE=F1, STATUS=’ OLD’ ) 

OPEN <2, FILE-F2, STATUS^’ NEW’ ) 

OPEN <3, FILE— F3, STATUS=’ NEW’ ) 

DO 1 1=1, N 

READ<1, *) < A < J, I ) , J= 1, M) 

IF <LQGY. NE- 1 ) GOTO 79 
DO 78 1 = 1, N 
DO 78 J=1,M 
A < J, I ) =AL0G1 0 < A < J, I) ) 

CONTINUE . 

NM=M 

CALL SVD<NM, M, N, A, SIGMA, . TRUE 


< A) ’ 
’ < A) 
< A) ’ 
’ < A) 
< A) ’ 
’ < A) 


’ ) ’ 

) F 1 
’ ) ’ 
)F2 
’ ) ’ 
)F3 
’ ) ’ 


ENTER INPUT FILNAM’ 


ENTER OUTPUT BASIC WAVEFORM FILNAM’ 


ENTER OUTPUT SVD COEFFICIENT FILNAM’ 


ENTER 1 FOR SORTING SV’ 


ENTER 1 FOR LOG Y’ 


U, . TRUE. , V, IERR, WORK) 


WRITE <*, ’ < A) ’ ) ’ VECTOR NORM OF U ARE:’ 

DO 501 J=l, N 

P=0. 

DO 502 1=1, M 
P=P+U<I, J>**£ 

WRITE <*, *) P 
CONTINUE 

WRITE < *, ’ < A) ’ ) ’ VECTOR NORM FOR V ARE:’ 

DO 503 J= 1 , N 

P=0. 

DO 504 1=1, N 
p=p+V<I, J>**£ 

WRITE <*, *) P 
CONTINUE 

WRITE <*, ’ <A) ’ ) ’ SINGULAR VALUES ARE:’ 

IF < IFST. EQ- I) THEN 
DO 777 1 = 1, N 
WK < I ) =S I GM A ( I ) 

CALL SORT <WK, I JK, N) 

ELSE 


DO 333 1 = 
I JK < I ) = I 


1 1 «* N 


D 

Line# 

1 


SO 



61 

C 


6£ 



63 



64 

£ 


65 


1 

66 


£ 

6^ 

3 


68 


1 

69 


2 

70 

666 


71 


1 

7£ 


£ 

73 


£ 

74 

201 


75 


1 

76 


£ 

77 

£0£ 


78 


1 

79 


1 

80 


1 

81 

4 


8£ 



83 



84 
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ENDIF 

WRITE (*, ’ (A) ’ ) ’ SORTING ORDER OF SINGULER VALUES:’ 

WRITE (*,•#) (SIGMA ( I JK ( I ) ) , 1 = 1, N) 

IF(IERR. NE.O) WRITE (*,£) I ERR 
FORMAT (’ TROUBLE. IERR=’ , 14) 

DO 3 1 = 1, .N 
DO 3 J=l, N 

CO (I, J)=V<I, J)*SIGMA(J> 

DO 666 1 = 1, N 
DO 666 J=1,N 
DD ( I, J) =CC ( J, I) 

DO £01 1 = 1, N 
DO £01 J=1,M 
A < J, I ) =U ( J, I JK ( I ) ) 

CONTINUE 

DO £02 1=1, N 

DO £0£ J=1,N 

CC ( I , J ) =DD < I JK < I ) , J ) 

DO 4 J= 1 , N 

WRITE (2, *) ( A ( I , J) , 1=1, M) 

WRITE (3,*) (CC ( I, J) , 1=1, N) 

CONTINUE 

WRITE (3,*) (SIGMA (IJK (I) ) , 1=1, N) 

STOP 

END 


Name 

Type 

Offset 

a 

AL0G10 

REAL 

39466 

CC 

REAL 

38498 

DD 

REAL 

38982 

FI 

CHAR*16 

57074 

F£ 

CHAR*16 

57090 

F3 

CHAR* 16 

57106 

I 

INTEGERS 

57130 

IERR 

INTEGERS 

57158 

IFST 

INTEGER*4 

57122 

IJK 

INTEGER*4 

2 

J 

INTEGER*4 

57138 

LOGY 

INTEGER*4 

57126 

M 

INTEGERS 

57066 

N 

INTEGER*4 

57070 

NM 

INTEGERS 

57154 

P 

REAL 

57166 

SIGMA 

REAL 

36898 

U 

REAL 

1698 

V 

REAL 

19298 

WK 

REAL 

50 

WORK 

REAL 

98 


Class 


INTRINSIC ‘ 



Q5 

SUBROUTINE SORT ( X , IJK, N) 


86 

DIMENSION XCl), IJK(l) 


87 

AM I N=X (1) 


88 

DO 1 1 = 1, N 

1 

89 

IF(X (I) . LT. AMIN) AMI N=X (I 

1 

90 1 

CONTINUE 


91 

AM I N=AM I N— 1 GO- 



- 
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92 

DO 2 1=1, N 



1 

93 

COLL FINDMAX (X, N, MAX) 



1 

94 

I JK < I ) =MAX 



1 

95 

X (MAX) =AMIN 



1 

96 2 

CONTINUE 




97 

RETURN 




96 

END 



Name Type 

Offset P Class 



AMIN REAL 

57254 



I 

INTEGER*^ 57258 



I JK 

INTEGER*^ 4 * 



MAX 

INTEGER** 57270 


» 

N 

INTEGERS 8 * 



X 

REAL 

o * 




99 C***************************** 




100 

SUBROUTINE FINDMAX (X, N, MAX) 




101 

REAL** X ( 1 ) 




102 

AMAX=X (1) 




103 

MAX=1 




104 

DO 1 1=2, N 



' 1 

105 

IF ( X < I ) . GT. AMAX ) THEN 



1 

106 

AMAX=X (I) 



1 

107 

MAX = I 



1 

108 s 

END IF 



1 - 

109 1 

CONTINUE 




110 

RETURN 




111 

END 



Name 

Type 

Offset P Class 



AMAX 

REAL 

57274 



I 

INTEGERS 57278 



MAX 

INTEGER*4 8 * 



N 

INTEGERS 4 * 



X 

REAL 

0 * 




112 

SUBROUTINE SVD (NM, M, N, A, W, MATU, U, MATV, V, IERR, 

RV1 ) 



113 

REAL*4 A (NM, N ) , W < N ) , U ( NM, N) , V (NM, N) , RV1 (N) 




114 

LOGICAL MATU, MATV 




115 

I ERR=0 




116 

DO lOO 1 = 1, M 



1 

117 

DO 100 J=1,N 



2 

118 

U(I, J)=A(I, J) 



2 

119 100 

CONTINUE 




120 

G=0. 




121 

SCALE=0. 




122 

AI\)0RM=0. 




123 

DO 300 1 = 1, N 



1 

124 

L=I + 1 



■» 

1 25 

RV1 <I)=SCALE*G 



i 

1 2S 

G=0. 



1 

127 

s=o 



1 

128 

5CALZ=0. 



i 

129 

Ir ( I. GT. M) GOTO 210 








2 

> 
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1 

130 


DO 120 K=I,M 

2 

131 

ro 

o 

SCALE=SCALE+ABS (U (K, I ) ) 

1 

132 


IF (SCALE- EQ. 0. 0) GOTO 210 

1 

133 


DO 130 K= I , M 

r» 

134 


U<K, I ) =*U <K, I ) /SCALE 

2 

135 


S=S+U<K, I)**2 

* 2 

136 

130 

CONTINUE 

1 

137 


F=U(I, I) 

1 

133 


G=-SIGN (SQRT <S) , F) 

' • i 

139 


H=F*G-S 

* l 

140 


U(I, I ) =F— G 


, 1 

141 


IF ( I - EQ. N) GOTO 190 


i 

142 


DO 150 J=L, N 


2 

143 


S=0. 


£ 

144 


DO 140 K=I,M 


145 

140 

S=S+U(K, I ) *U <K, J) 

[ 2 

146 


F=S/H 

2 

147 


DO 150 K= I , M 

vi 

143 


U (K, J) =U <K, J) +F*U <K, I) 

3 

149 

150 

CONTINUE 

1 

150 

190 

DO 200 K=I,M 

; 2 

151 

200 

U (K, I>=SCALE*U(K, I) 

> 1 

152 

210 

W < I ) =SCALE*G 

1 

153 


G=0. 

i 1 

154 


S=0. 

1 

155 


SCALE=Q. 

] 1 

156 


IF(I. GT. M. OR. I. EQ. N)GOTO 290 

I 1 

157 


DO 220 K=L, N 

2 

158 

220 

SCALE=SCALE+ABS(U(I, K> > 

’ 1 

159 


IF (SCALE- EQ- 0. 0) GOTO 290 

; 1 

160 


DO 230 K=L, N | 

2 

161 


U(I,K)=U(I,K) /SCALE 

2 

162 


S=S+U(I, K>**2 

2 

163 

230 

CONTINUE 

1 

164 


F=U<I* L) 

1 

165 


G=— SIGN (SQRT <S) , F) 

1 

166 


H=F*G-S 

1 

* 167 


U(I, L)=F-G 

’ 1 

168 


DO 240 K=L,N 

2 

169 

240 

RV1 (K) =U< I, K) /H 

1 

170 


IF (I. EQ. M) GOTO 270 

1 

171 


DO 260 J=L, M 

1 o 

172 


S=0. 

2 

173 


DO 250 K=L,N 

3 

174 

250 

S=S+U(J, K)*U(I, K) 

p 

175 


DO 260 K=L,N 


176 


U ( J, K) =U ( J, K) +S*RV1 (K) 

3 

177 

260 

CONTINUE 

1 

178 

270 

DO 280 K=L, N 

a 

179 

280 

U ( I , K) =SCALE*U ( I , K) 

1 

180 

290 

ANORM=AMAX 1 ( ANORM, ABS < W ( I ) >+ABS(RVl (I) ) ) 

1 

181 

300 

CONTINUE 


182 


IF (. NOT. NATO) GOTO 410 


183 


DO 400 11=1, N 

1 

104 


I =N+ 1 - I I 

1 

185 


IF(I.EQ.N)GOTQ 390 

-• 

186 


IF ( G . c.0 - 0- O) GOTO 360 

X 

1 87 


DO 320 J=L, N 

2 

188 

320 

V(J, I) = (U(I, J)/U<I,L) )/G 


1 
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1 

189 


DO 350 J=L, N 


2 

190 


S=0. 


2 

191 


DO 340 K=L, N 


3 

i9a 

340 

S=S+U(I, K)*V<K, J) 


2 

193 


DO 350 K=L, N 


3 

194 


V<K, J)=V(K, J)+S*V<K, I) 


3 

195 

350 

CONTINUE 


1 

196 

360 

DO 380 J=L,N 


2 

197 


V<I, J)=0. 


2 

198 


V(J,I)=0. 


2 

199 

380 

CONTINUE 

* 

1 

200 

390 

V(I,I)-1. 


1 

201 


G=RV 1 (I) 


1 

202 


L=I 


a 

£03 

400 

CONTINUE 



204 

410 

IF < . NOT. MATU) GOTO 510 



£05 


MN=N 



206 


IF <M. LT. N>MN=M 



207 


DO 500 11=1, MN 


1* 

208 


I=MN+1-I I 


1 

209 


L-I + l 


1 

21C 


G=W ( I ) 


1 

21 1 


IF (I. EG?. N) GOTO 430 


1 

212 


DO 420 J=L, N 


£ 

£13 

4£0 

U(I, J) =0. 


1 

214 

430 

IF ( G. EQ. 0. ) GOTO 475 


1 

£15 


IF < I . EQ. MN) GOTO 460 


1 

£16 


DO 450 J=L, N 


2 

£17 


S=0. 


2 

£18 


DO 440 K=L,M 


3 

219 

440 

S=S+U(K, I ) *U (K, J) 


2 

220 


F= (S/U < I , I) ) /G 


2 

221 


DO 450 K=I,M 


3 

£22 


U<K, J>=U<K, J)+F*U(K, I) 


3 

£23 

450 

CONTINUE 


1 

£24 

460 

DO 470 J=I,M 


2 

225 

470 

U(J, I)=U(J, I) /G 


1 

£26 


GOTO 490 


1 

£27 

475 

DO 480 J=I , M 


£ 

228 

480 

U(J, I ) =0. 


1* 

£29 

490 

U<I, I)=U(I, I)+l. 


1 

£30 

500 

CONTINUE 



£31 

510 

DO 700 KK=1,N 


1 

£32 


K 1 =N— KK 


1 

233 


K=K1+1 


1 

£34 


I TS=0 


1 

235 

520 

DO 530 LL= 1 , K 


£ 

£36 


L 1 =K-LL 


£ 

237 


L=L1 + 1 


2 

238 


IF (ABS (RV1 (L) ) +ANORM. ED. ANORM) GOTO 565 


2 

239 


IF ( ABS ( W <L 1 ) ) +ANQRM. EQ. ANORM ) GOTO 540 


2 

240 

530 

CONTINUE 


1 

241 

540 

c=o 


-4 

X 

242 


S= 1 . 


l 

243 


DO 560 I=L,K 


iZ. 

244 

- 

F=3*RV1 ( I ) 


2 

245 


RV1 ( I ) =C*RV1 < I ) 


2 

246 


IF (ABS (F) + ANORM. EQ. ANORM) GOTO 565 


CL 

247 


G=Ui ( I ) 








. 

i 



• 



? 
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i 

2 

248 

H=SQRT<F*F+G*G) 



* 

2 

249 

W ( I ) =H 




2 

250 

C=G/H 




2 

251 

S=— F/H 




2 

252 

I F ( . NOT. MATU ) GOTO 560 




2 

253 

DO 550 J=?1,M 



I 

3 

254 

Y=U ( J, LI) 



V 

3 

255 

Z=U(J, I) 



i 

3 

256 

U(J, L1)=Y*C+Z*S 



j 

3 

257 

U<J, I)=-Y*S+Z*C 



4 

3 

258 550 

CONTINUE 


' i 

V 

\ 

2 

259 560 

CONTINUE 


( 

> 

3 

1 

260 565 

Z-W (K) 


! 

• i 

1 

361 

IF (L. EQ. K) GOTO 650 


t 


1 

262 

IF ( ITS. EQ. 30) GOTO 1000 




1 

263 

ITS=ITS+1 




1 

264 

X=W(L) 




1 

265 

Y=W(K1) . 



f 

1 

266 

G=RV 1 (Kl) 



i 

i 

1 

267 

H=RV1 <K) 



"1 

1 

268 

F=< < Y-Z) * < Y+Z ) + (G— H) * <G+H) )/<2. *H*Y) 



* 

1 

269 

G—SQRT <F*F+1. ) 



I 

! 

1 

270 

F=< <X-Z)* <X + Z) +H*(Y/ (F+SIGN <G f F) >-H) ) /X 




1 

271 

Ol. 



J 

1 

£7£ 

S=l. 



t 

1 

273 

DO 600 I 1=L, Kl 



5 

$ 

2 

£74 

1 = 114-1 



; 

2 

275 

G=RV 1(1) 



i 

2 

276 

Y=W < I ) 



{ 

2 

277 

H=S*G 




2 

278 

G=C*G 



t 

2 

279 

Z=SQRT<F*F+H*H) 



j 

2 

280 

RV1 <ID=Z 



i 

2 

281 

C=F/Z 



i 

2 

282 

S=H/Z 




2 

283 

F=X*C+G*S 



- 

2 

284 

G=-X*5+G*C 



* 

2 

285 - 

H=Y*S 



< 

o 

c 

286 

Y=Y*C 




2 

287 

IF < . NOT. MATV) GOTO 575 




2 

288 

DO 570 J= 1 , N 



j 

3 

289 

X=V< J, ID 




3 

290 

Z=V(J, I) 




3 

291 

V(J, I1)=X*C+Z*S 




3 

292 

V(J, I)=-X*S+Z*C 



r 

3 

293 570 

CONTINUE 



j 


294 575 

Z=SQRT <F*F+H*H) 




2 

295 

W(I1)=Z 



* < 

2 

296 

IF(Z. EQ. 0. 0)G0T0 580 




2 

297 

C=F/Z 




2 

298 

S=H/Z 




2 

299 580 

F=C*G+S*Y 




2 

300 

X=-S*G+C*Y 




2 

301 

IF (. NOT. MATU) GOTO 600 




2 

302 

DO 590 J=1,M 




«J 

303 

Y=U ( J, ID 





304 

Z=U(J, I) 




3 

305 

U(J, 1 1 ) =Y*C+-Z*S 




o 

306 

U(J, I ) =— Y*S+Z*C 





- - ■ 

- 


- 


d: 

Line# 

i 

7 

3 

o07 

590 

CONTINUE 

c. 

308 

600 

CONTINUE 

X 

309 


RV1 (L) =0. 

1 

310 


RV1 (K) =F 

1 

31 1 


W(K)=X 

1 

312 


GOTO 520 

1 

313 

650 

IF ( Z. GE. 0.0) GOTO 700 

1 

314 


W (K) =— Z 

1* 

315 


IF (. NOT. MftTV) GOTO 700 

L. 

316 


DO 690 J= 1 , N 

2‘ 

317 

690 

V < J, K ) =— V < J, K) 

1' 

318 

700 

CONTINUE 


319 
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320 

1000 
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Name 

Type 

ft 

REAL 

ftBS 

ftMftXl 

ftNORM 

REAL 

C 

REAL 

F' 

REAL 

G- 

REAL 

H 

REAL 

i: 

INTEGER*4 

n 

INTEGER*4 

IERR 

INTEGER*^ 

1 1 

INTEGER*^ 

ITS 

INTEGERS 

J 

INTEGER*^ 

K. 

INTEGER*^ 

K 1 

IIMTEGER*4 

KK 

INTEGER*^ 

L- 

INTEGER*^ 

LI 

INTEGER*^ 

LL 

INTEG£R*4 

M 

INTEGER*4 

MftTU 

LGGICAL*4 

MftTV 

L0GICAL*4 

MN 

INTEGERS 

N 

INTEGERS 

NM 

INTEGER*4 

RV1 

REAL 

S 

REAL 

SCftLE 

REAL 

SIGN 

SORT 

U 

REAL 

V 

REAL 

w 

REAL 

X 

REAl 

Y 

REAl 

z 

REAl. 


Offset P Class 
12 * 

INTRINSIC 

INTRINSIC 

57310 
57478 
57 o3Q 
57302 
57342 
57286 . 

57502 
36 * 

57390 

57462 

57294 

57326 

57458 

57450 

57318 

57474 

57466 

4 * 

20 * 

28 * 

57418 

8 * 

0 * 

40 * 

57322 

57306 

INTRINSIC 

INTRINSIC 

24 * 

32 
16 
57498 
57490 
5749 ** 
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FLOW DYNAMIC ENVIRONMENT 
DATA BASE DEVELOPMENT FOR THE SSME 

C. V. Sundaram 

Wyle Laboratories, Huntsville 
Abstract 


This paper describes the studies being carried 
out on the fluid flow-induced vibration of the Space 
Shuttle main engine (SSME) components. This study is 
being carried out with a view to correlating the fre- 
quency characteristics of the pressure fluctuations in 
a rocket engine to its operating conditions and geom- 
etry. An overview of the data base development for 
SSME test firing results and the interactive computer 
software used to access, retrieve, and plot or print 
the results selectively for given thrust levels, 
engine numbers, etc., is presented. The various sta- 
tistical methods available in the computer code for 
data analysis are discussed. Plots of test data, 
nondimensionalized using parameters such as fluid flow 
velocities, densities, and pressures, are presented. 
Comparative studies of these results with results 
available in the literature are made. Correlations 
between the resonant peaks observed at higher frequen- 
cies in power spectral density plots with pump geome- 
try and operating conditions are discussed. Finally, 
an overview of the status of the investigation is pre- 
sented and future directions are discussed. 

Introduction 


The development efforts expended on the main 
engine of the Space Shuttle has resulted in the evolu- 
tion of highly complex machinery which is capable of 
generating high thrust levels. Even from the earliest 
days of spacecraft engine development, achieving the 
highest possible specific thrust, based on engine 
weight or overall size, has been identified as one of 
the prime goals. The earlier spacecraft, however, 
were equipped with power sources which were designed 
for just one mission. The advent of the concept of 
reusable space vehicles added a new dimension to the 
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engine design, that they not only satisfy the require- 
ments on thrust but also that they be capable of 
carrying out multimissions into outer space with, 
ideally, the least amount of maintenance and component 
replacements between missions. This has imposed a 
greater emphasis on preventing failures and prolonging 
the fatigue life expectancy of the SSME during opera- 
tion. Furthermore, specific knowledge relating to the 
frequency characteristics of the pressure fluctua- 
tions in rocket engine systems is not well understood. 
It is important to understand the functional relation- 
ships between pressure fluctuations in these systems 
in terms of the geometrical and performance charac- 
teristics because of their effect on hardware 
lifetime. 

The SSME has been subjected to extensive hot 
firing and flow tests. During these tests, system 
failures and malfunctions have occurred from time to 
time. These failures range from wear on component 
bearings and flow- induced failures to explosions due 
to intense pressure fluctuations and accompanying 
dynamic stresses in the pumps, valves, and/or propel- 
lant lines. The large volume of data on pressure 
fluctuations, strain, and acceleration obtained from 
these tests provide a good basis for studies on the 
flow- induced vibrations taking place in the SSME. 

The sheer volume of test data available makes it 
absolutely essential that a computerized data base 
management system be employed to perform any meaning- 
ful analysis. Wyle Laboratories is engaged in the 
development of such a data base management software. 
This paper provides an overview of the work being per- 
formed on the software development and discusses the 
method of approach suggested for analyzing the data. 

. SSME Data Base Management 

The SSME data base management and analysis system 
developed at Wyle Laboratories has been designed for 
use on an Interdata 8/32 computer system and other 
computer systems utilizing similar file name 
structures. 


Data Filing System 


The data base management system is designed to 
accommodate two types o£ data: time histories and 
power spectral density or other frequency domain data. 
The time history (rms pressure or thrust) data files 
consist of a pair of values in each record, the first 
value being the number of seconds into the test and 
the second being the magnitude. The file may consist 
of any number of data points (records) , but most of 
the routines presently have an upper limit of 500 
points. 

Each PSD file consists of a single real number in 
each record. These are the spectral magnitudes. 
Since the PSD test data is recorded using constant 
bandwidth filters, the frequencies corresponding to 
each PSD value is automatically calculated by the pro- 
gram when the frequency range is entered. A consider- 
able saving in computer disc storage space has been 
realized by this method because of the large amount of 
PSD data involved. 

Data entry into files can be performed through a 
digitizing tablet. When the volume of data is large, 
however, use of a digitizing tablet can be very time 
consuming. Bence the data in this case is directly 
read from a digital tape using the dedicated magnetic 
tape drive units. File names are automatically 
assigned by the software based on test numbers, engine 
numbers, channel descriptors, etc., to uniquely iden- 
tify each data set. Figure 1 provides an insight into 
the operational sequences and capabilities of the SSME 
data base management and analysis software. 

Figure 2 shows a- plot of power versus time, and 
figure 3 shows the rms pressure variations with time 
for the low pressure fuel pump. These plots were 
generated using the data retrieval system described 
above. Figures 4 shows a PSD plot generated from the 
digitized data for the high pressure oxidizer pump. 
The input to the computer to generate these plots are 
the file names for the time history plots and the test 


number/ time slice/ and maximum frequency for the PSD 
plots. 


Statistical parameters/ such as mean f standard 
deviation/ rms values/ third and fourth central 
moments/ and minimum and maximum values and range/ can 
be calculated and plotted with ease for any number of 
data sets selected based on engine/ HPFP, and HPOP 
numbers; thrust level; and channel descriptors. The 
abscissa and the ordinate can have any combination of 
linear and logarithmic scales/ and any statistical 
parameter can be plotted against any other or a fre- 
quency scale. The uniqueness of the set of PSDs 
selected for statistical analysis can be controlled by 
specifying or leaving out/ in any combination/ the 
engine/ HPOP/ and HPFP numbers. 

A detailed description of all the capabilities of 
the data base management and analysis software is pre- 
sented in reference 1. 

General Character of Turbomachinery Noise 
Background 

The jet noise theory gained importance with 
increasing use of jet engines in aircraft. It is well 
known that the acoustic power generated by a jet 
exhaust is strongly dependent* on the exit velocity 
(Ligh thill's 8th power law z ) . However/ during 
approach configuration of the aircraft/ when the jet 
noise level is generally low, the engine compressor 
noise was observed to be a major noise source. The 
study of compressor noise/ especially the multistage 
axial flow types/ has acquired great importance during 
the last 25 years. 

Noise Sources in Axial Flow Compressors 

In general/ the noise spectra of different types 
of compressors exhibit similar gross characteristics. 
In this paper/ we limit our attention to axial flow 
compressors because it is believed that their noise 
characteristics should be more relevant to the noise 


generation mechanism of the axial flow pumps on the 
SSME. 


The general form of an acoustic spectrum of an 
axial flow compressor would indicate broadband noise 
extending over a wide range of frequencies. Super- 
imposed on this are a number of discrete peaks which 
represent the blade passage frequencies and their 
harmonics. 

The broadband noise component has been attributed 
to various mech^pisms. At least three that merit con- 
sideration are 

1. The random force fluctuations due to a moving 
flow on the surface. These fluctuations can be inter- 
preted to act as acoustic sources located within the 
turbulent boundary layer. 

2. The shedding of vortices from the trailing edges 
of a body immersed in a moving fluid. This vortex 
shedding imparts to the body fluctuating lift forces 
that are periodic at low Reynolds number flows and 
became randan as the Reynolds number increases. (It 
is believed that the periodic component of the shed- 
ding does not disappear completely.) This random exci- 
tation force of vortex shedding is also considered to 
be instrumental in the broadband characteristic of the 
acoustic spectrum. 

3. A turbulent flow at the inlet to the compressor 
sets up fluctuating forces due to flow incidence on 
the surface. The perturbation velocity component in a 
turbulent flow is random in nature/ and hence the 
acoustic energy due to the fluctuating forces can be 
expected to have a broadband nature. 

Another significant broadband component in the 
case of the SSME would be the structurally transmitted 
secondary vibrations due to the close proximity of the 
components/ all having highly turbulent flow taking 
place within them. A challenging task in the analysis 
of SSME component vibration is the development of a 
prediction scheme that will provide quantitative 
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estimates of the acoustic energy due to flow- induced 
and structurally transmitted vibrations in a frequency 
spectrum at various locations on the SSME. 

The pure tones, which are a result of the spin- 
ning modes of the pressure field, exhibit decaying 
trends downstream of the rotor when the impeller tip 
Mach number is below critical, as in the case of a 
subsonic rotor. For a rotor-stator combination, how- 
ever, a number of spinning modes are generated, many 
of which spin at Mach numbers above the critical tip 
Mach numbers, and thes| modes exhibit strong propaga- 
tion characteristics. 

SSME Flow- Induced Vibration Characteristics 


In this section, preliminary studies conducted on 
pressure spectra obtained from SSME firing tests are 
discussed. Fluid flow variables, which may be of 
lesser importance in conventional axial flow compres- 
sors, but which may be of greater significance in the 
SSME flow environment, are highlighted. 

The initial efforts in the data analysis have 
been directed toward obtaining nondimensional PSD as a 
function of Strouhal number. The power spectral den- 
sities, which have dimensions of mean square pressure 
per hertz, were normalized as follows: 

% m V f)/P f ^ ° (1) 

where S is nondimensional PSD, 0(f) is the measured 
PSD, is the local fluid density, V is the mean 
velocity, and D is a characteristic length. Strouhal 
number is defined as 

S fc - fD/V (2) 
where f is the frequency. 

Figures 5 and 6 show nondimensionalized PSD plot- 
ted as a function of Strouhal number for the high 
pressure oxidizer pump discharge pressure. Similar 
analyses have been carried out to estimate the 
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pressure fluctuations due to turbulent flow of water 
at the wall of a duct. Figure 7, reproduced from 
reference 5, shows the pressure fluctuation dependence 
on Strouhal number. Figures 5, 6, and 7 indicate 
similar trends in the variation of nondimensionalized 
PSD for increasing Strouhal number. The scaling 
applied to the SSME data is based on the assumption 
that the pressure fluctuations are flow induced. This 
assumption is somewhat oversimplified as it does not 
address the various other noise-generating mechanisms 
present during turbomachinery operation. 

6 

Dimensional analysis studies of blower noise 
provide a functional relationship, 

S - F(D,N,Q,fi,p,fi, C) (3) 

where E is the acoustic power/hertz; D, the blower 
diameter; N, the rpm; Q, the flow rate; f, the fre- 
quency; p, the fluid density; p, the viscosity; and C, 
the acoustic speed. A similar analysis on the SSME 
turbopump would require accounting for the acoustic 
propagation characteristic of the structure. 

Another important factor that should be con- 
sidered is the fluid properties. The fluid flow 
conditions in the SSME undergo extreme variations in 
temperature and pressure. The mechanism of propaga- 
tion of disturbances in fluids at extremely low tem- 
peratures and very high pressures requires further 
studies for proper mathematical modeling. Since the 
fluid flow conditions experience wide variations in 
the SSME, it is reasonable to expect variations in 
local propagation velocities and hence local Mach 
numbers. 

The studies on the flow dynamic environment of 
the SSME being carried out at wyle Laboratories are 
directed toward applying the existing theories on 
noise generation mechanism in turbomachinery after 
modifications to reflect the flow environment present 
in the SSME. An approach that would involve develop- 
ing an empirical formulation rather than a very rigor- 
ous and purely analytical solutions seems to indicate 


promise. The basic idea involves modeling the noise 
generation mechanism due to turbulent flow and the 
pressure pulses generated by the spinning modes o£ the 
pressure field due to rotor-stator interactions sepa- 
rately and combining the two components to obtain the 
total frequency spectrum. That an empirical approach, 
though somewhat complicated, is not totally impossible 
has been amply demonstrated through similar studies 
performed by Van Niekerk and others. 
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Figure 1. SSME data base management and analysis 
software operational sequence and capabilities 
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Figure 2. SSME power level variations as a function 
of time 



Figure 3. rras pressure time history at the inlet of 
the low pressure fuel pump 
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Figure 4. Power spectral density of the high pressure 
oxidizer pump discharge pressure 5 seconds after SSME 
startup 
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Figure 5. Nondimensionalized frequency spectrum for 
high pressure oxidizer pump discharge pressure, 5 
seconds after SSME startup 
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Figure 6. Nondimensionalized frequency spectrum for 
high pressure oxidizer pump discharge pressure 400 
seconds after SSME startup 



Figure 7. Frequency spectrum of turbulent wall pres- 
sure field. (O), 267 in/sec; (A) 450 in/sec; (□) 520 
in/sec 3 
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